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(i) 


ABSTRACT 


The  relativistic  gravitational  collapse  of  an  object  is  usually 
considered  in  four  stages:  instability,  implosion,  horizon,  and  singular¬ 
ity.  In  addition,  the  collapse  itself  falls  into  one  of  three  categories: 
spherical,  nearly-spherical ,  or  highly  nonspherical .  This  thesis  is 
devoted  to  the  presentation  of  three  studies,  each  of  which  deals  with  a 
particular  phase  of  this  complex  problem. 

In  the  first  study  we  present  some  fairly  elementary  results, 
related  to  the  instability  stage  of  spherical  collapse.  For  simplicity 
we  deal  with  spherical  shells.  A  solution  of  the  Einstein  equations  is 
derived,  representing  a  thin  spherical  shell  of  charged  fluid,  falling  in 
a  spherically  symmetric  field  due  to  mass  and  charge  at  its  centre.  No 
restrictions  are  placed  on  the  equation  of  state.  We  integrate  the  equa¬ 
tions  of  motion  to  find  the  law  of  conservation  of  total  energy,  and  we 
use  it  to  study  the  equilibrium  states  of  the  system,  and  their  stability 
against  collapse.  We  find  under  reasonable  assumptions,  that,  given  the 
entropy  and  the  equation  of  state,  there  is  a  maximum  equilibrium  mass,  and, 
corresponding  to  it,  a  critical  radius,  inside  of  which  instability  sets  in. 
For  uncharged  bodies,  these  results  completely  parallel,  and  serve  as  a 
simple  illustration  of,  the  much  more  complicated  analyses  needed  for  fluid 
spheres. 


The  second  study  is  concerned  with  the  implosion  and  horizon 


stages  of  nearly-spherical  collapse.  Two  idealized  collapse  models  (again 
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(ii) 


thin  shells),  involving  a  scalar  monopole  and  a  magnetic  dipole,  are 
considered,  treating  departures  from  sphericity  as  small  perturbations. 
Radiative  leakage  (largely  downwards  through  the  Schwarzschild  horizon) 
causes  the  externally  observable  asymmetries  to  decay  to  zero  in  an 
oscillatory  fashion.  These  results  have  significant  consequences  for 
astrophysics;  they  imply  in  particular  that  a  "black  hole"  cannot  be  a 
source  of  synchrotron  radiation. 

In  the  third  study  we  deal  mainly  with  the  horizon  stage  of 
highly  nonspherical  collapse,  and  we  consider  static  scalar  fields.  We 
prove  the  following  theorem.  Every  zero-mass  scalar  field  which  is  gravi¬ 
tationally  coupled,  static  and  asymptotically  flat,  becomes  singular  at  a 
simply-connected  event  horizon.  In  the  special  case  where  the  gravitational 
coupling  of  the  scalar  energy  density  is  neglected,  the  solutions  are 
computed  explicitely.  Some  properties  of  the  singular  event  horizons  are 
also  discussed. 
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CHAPTER  I 

Introduction 


§1.1  Spherical  Collapse 

The  story  of  relativistic  gravitational  collapse  really  began 
in  1916,  when  the  normally  formidible  Einstein  equations  were  shown,  by 
Schwarzschild  [1],  to  yield,  in  the  static  spherically  symmetric  vacuum 
case,  a  beautifully  simple  solution, 

ds  =  -(1 — )dt  +  (1 — )  dr  +  r  (d0  +  sin  0  d<j>  )  .  (1.1) 

(We  choose  the  signature  to  be  -  +  +  +  and  our  units  such  that 
G  =  c  =  1  .) 

Several  of  the  well-known  properties  of  the  metric  form  (1.1) 
are  worth  noting. 

(i)  Recall  that  (1.1)  is  static.  It  can  be  shown  that  "any 

spherically  symmetric  vacuum  field  is  static".  This  result 
is  known  as  Birkhoff's  Theorem  [2],  and  in  particular  it 
tells  us  that,  even  in  the  case  of  a  star  undergoing  catas¬ 
trophic  spherical  collapse,  the  external  line  element  is 
just  the  Schwarzschild  metric  (1.1). 

(ii)  The  apparent  singularity  at  r  =  2m  is  not  a  physical 

singularity,  but  merely  a  pathology  of  the  particular  coordi¬ 
nate  system  chosen  -  in  fact,  by  transforming  to  different 
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coordinates,  or  examining  scalars  formed  from  the  curva¬ 
ture  tensor,  it  is  easily  seen  that  the  space-time  is 
perfectly  well-behaved  at  r  =  2m  [3].  On  the  other  hand, 
the  r  =  2m  hypersurface  does  have  special  properties.  It 
is  an  "event  horizon"  [4],  i.e.  it  is  the  inner  boundary  of 
the  set  of  events  which  can  be  connected  to  spatial  infinity 

by  future  timelike  lines.  Also,  every  spacelike  2-sphere 
2 

S  (r  =  const.  <  2m  ,  t  =  const.)  inside  the  r  =  2m 
hypersurface  is  "trapped",  i.e.  both  systems  of  null 

2 

geodesics  (outgoing  as  well  as  incoming) ,  which  meet  S 
orthogonally,  converge  locally  in  future  directions  [5], 

(iii)  The  singularity  at  r  =  0  is  a  true  physical  one  [6],  i.e. 

a  free-falling  observer  would  be  torn  apart  by  tidal  gravi¬ 
tational  forces  near  r  =  0  . 

Since  the  gravitational  surface  r  =  2m  is  not  physically 
singular,  there  is  no  physical  reason  why  a  star  should  not  fall  inside 
it.  At  the  same  time,  once  a  star  is  inside  its  own  gravitational  radius, 
it  can  never  reemerge  (property  (ii)),  hence  collapse  is  irreversible. 

In  this  sense  then,  general  relativity  predicts  the  inevitability  of 
spherical  stellar  collapse,  and  these  predictions  have  been  verified,  for 
certain  idealized  models,  by  several  calculations  [7]. 

The  overall  phenomenon  of  spherical  gravitational  collapse  has 
been  portrayed  by  Thorne  [8]  as  having  four  key  stages: 
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(i)  Instability.  At  a  certain  stage  in  its  evolution,  the 

star  runs  out  of  nuclear  fuel,  which  affects  the  equation 
of  state  and  thus  leads  to  instability. 

(ii)  Implosion.  In  a  very  short  time  the  instability  causes 

collapse  to  begin,  with  the  dense  core  falling  more  rapidly 
and  leaving  an  outer  envelope  to  trail  behind. 

(iii)  Horizon.  The  stellar  surface  reaches  its  gravitational 

radius  in  a  finite  time,  as  measured  by  a  comoving  observer. 
However,  to  a  distant  observer,  the  star  takes  infinite  time 
to  approach  r  =  2m  ,  i.e.  this  surface  is,  for  him,  an 
event  horizon.  This  horizon  becomes  the  boundary  of  a  "black 
hole",  a  region  of  space-time  that  can  never  communicate  with 
the  "outside  world". 

(iv)  Singularity .  The  star  continues  to  collapse  and  soon  reaches 
the  singularity  at  r  =  0  ,  corresponding  to  infinite  density 
and  infinite  tidal  gravitational  forces. 

§1.2  Near ly-Spheri cal  Collapse. 

A  fundamental  problem  still  remains.  Does  our  present  descrip¬ 
tion  of  idealized  spherical  collapse  also  fit,  to  any  degree  of  accuracy, 
the  case  of  a  more  realistic,  non-spherically  collapsing  star?  In  other 
words,  will  a  star,  whose  collapse  is  nonspherical  because  of  rotation, 
magnetic  fields,  asymmetric  density,  etc.,  behave  qualitatively  like  the 
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idealized  spherical  case?  Or  is  it  possible  that  these  complications 
might  alter  the  spherical  collapse  picture  by  stopping  or  reversing  the 
collapse  [9]  or  by  eliminating  the  event  horizon  [10],  thus  leaving  the 
r  =  0  singularity  visible  to  all? 

A  complete  solution  to  this  problem  still  is  not  known, 
although  several  recent  results  have  provided  us  with  a  good  picture  of 
collapse,  in  the  special  case  of  small  perturbations  from  sphericity.  We 
will  briefly  outline  these  developments  now. 

In  1963,  we  were  again  very  fortunate  when  the  Einstein  equa¬ 
tions  miraculously  yielded  another  simple  solution,  namely,  that  of  Kerr 
[11],  for  stationary  axi-symmetric  vacuum  fields.  This  was  generalized 
to  include  the  case  of  charge,  a  short  time  later  [12]. 

Early  studies  of  perturbed  spherical  collapse,  while  lacking 
detail,  indicated  that:  (i)  the  external  field  of  a  collapsing  star 
should  become  asymptotically  stationary  [13]; 

(ii)  if  we  assume  (i) ,  then  in  order  to  avoid 
singularities  [14],  all  externally-viewed  asymmetries  must  somehow 
disappear  [15]  as  the  collapsing  object  approaches  the  horizon. 

More  recently,  these  results  have  been  largely  confirmed. 

Novikov  [16]  proves  that,  in  fact,  the  perturbations  do  not  become  singular 
near  the  horizon.  De  la  Cruz,  Chase,  and  Israel  [17]  numerically  trace  the 
perturbed  collapse  of  several  simple  idealized  models,  and  show  that,  as 
predicted,  the  perturbations  do  not  blow  up,  but  are  "radiated  away"  as 
gravitational,  electromagnetic,  or  scalar  waves  -  some  going  off  to 
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infinity,  and  the  rest  falling  in  through  the  horizon. 

The  definitive  contribution  in  this  area  is  undoubtedly  the 
beautiful  work  of  Price  [18].  He  also  carries  out  numerical  analyses 
of  perturbed  collapse,  and  is  the  first  to  succeed  in  showing  how  the 
"tails"  of  the  radiated  waves  die  out,  for  large  t  .  Even  more  impor¬ 
tant,  he  proves  a  theorem  that  explains  exactly  why  the  perturbations 
are  radiated  away. 

We  can  therefore  summarize  nonspherical  gravitational  collapse, 
in  the  case  of  small  perturbations,  in  the  following  way: 

(i)  Nearly-sherical  collapse  does  behave  qualitatively  like 
idealized  spherical  collapse  (except  possibly  for  the 
"singularity"  stage,  about  which  little  is  known).  In 
other  words,  perturbed  spherical  collapse  also  fits  into 
the  standard  mold  -  instability,  implosion,  horizon. 

(ii)  Price's  Theorem:  "In  relativistic  gravitational  collapse 
with  small  nonspherical  deformations,  anything  that  can  be 
radiated,  will  be  radiated  away  completely."  This  theorem 
holds  for  zero-rest-mass  fields  of  arbitrary  integer  spin. 

It  is  well  known  that  all  multipole  moments  of  order  t  >_  s 
(for  a  massless  spin-s  field)  can  be  radiated  away  -  the 
rest  cannot,  i.e.  they  are  "conserved".  According  to  Price's 

i 

theorem,  then,  the  end  result  of  nearly-spherical  collapse 
must  be  a  "black  hole"  exhibiting  no  scalar  field,  only  a 
monopole  electromagnetic  field  (charge),  only  monopole  and 
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dipole  gravitational  fields  [19],  and  so  on. 

§1.3  Highly  Nonspherical  Collapse  [20] 

The  question  of  whether  or  not  black  holes  still  are  formed, 
when  the  collapse  is  of  a  highly  nonspherical  nature,  has  become  a 
widely  investigated  problem  in  recent  times  -  a  problem  not  yet  fully 
resolved.  Nevertheless,  mounting  evidence,  while  still  not  conclusive, 
points  to  the  following  results: 

(i)  "The  Black  Hole  Conjecture".  A  black  hole  forms  when  and 

only  when  a  mass  m  is  somehow  compressed  until  its  circum¬ 
ference  in  every  direction  is  C  <  4Trm  . 

(ii)  "A  black  hole  has  no  hair".  This  colourful  expression  of 
J.A.  Wheeler  memorably  paraphrases  the  conjecture  that, 
after  collapse,  the  horizon  will  eventually  settle  down  to 
a  stationary  state,  whose  external  field  is  given  by  the 
"charged  Kerr"  metric.  In  other  words,  a  black  hole,  and 
hence  its  exterior  geometry,  is  ultimately  determined  by 
just  three  "conserved  quantities"  -  the  mass  m  ,  charge  e  , 
and  angular  momentum  J  that  fall  into  the  hole.  Other 
factors,  such  as  the  collapsing  body's  asymmetries  (e.g. 
gravitational  quadrupole) ,  or  locally  conserved  quantities 
(e.g.  baryon  number),  have  no  net  effect  on  the  external 
field  [21]. 
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We  will  briefly  review  the  theorems  and  analyses  that  tend 
to  support  these  conjectures.  The  main  evidence  in  support  of  (i)  is 
naturally  its  well-established  truth  for  spherical  and  almos t-spherical 
collapse.  Several  examples  of  highly  nonspherical  collapse  have  been 
demonstrated,  in  which  no  horizon  is  formed  [22],  However,  these 
examples  are  pathological  in  nature,  and  the  absence  of  a  horizon  seems 
due  to  the  fact  that  the  collapse  does  not  occur  in  all  directions 
(hence  the  inclusion  of  the  phrase  "in  every  direction"  in  the  black 
hole  conjecture) . 

The  plausibility  of  conjecture  (ii)  is  enhanced  by  a  series 
of  rigourously  established  results,  concerning  static  and  stationary 
event  horizons.  We  can  summarize  the  static  (J  =  0)  results  as 
follows : 

(i)  The  Schwarzschild  fields  with  m  >_  0  are  the  only  static, 
vacuum  (e  =  0)  solutions  which  are  asymptotically  flat 
and  possess  simply-connected  equipotential  2-surfaces 
(gQQ  =  const.)  ,  and  a  nonsingular  horizon  gQQ  =  0  [23]. 

(ii)  The  Reissner-Nordstrom  solutions  with  m  >_  |e|  occupy  the 
same  unique  position  in  the  class  of  all  static  electrovac 
solutions  [24] . 

(iii)  Every  zero-mass  scalar  field  which  is  gravitationally 

coupled,  static  and  asymptotically  flat,  becomes  singular 
at  a  simply-connected  event  horizon  [25] . 
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From  these  theorems  we  would  like  to  conclude  that  "a  nonro¬ 


tating  black  hole  has  no  hair",  i.e.  it  is  completely  determined  by  the 
mass  m  and  charge  e  (if  any)  that  goes  down  the  hole.  Unfortunately, 
we  cannot  draw  this  conclusion,  for  two  reasons.  First,  the  topological 
requirement  that  the  equipotential  2-surfaces  be  simply-connected,  is 
overly  restrictive,  and  should  be  dropped.  Second,  the  demand  for  static 
solutions  is  too  strong  and  must  be  replaced  by  the  weaker  conditions 
"stationary  plus  J  =  0  ".  A  definite  improvement,  in  the  uncharged  case, 
has  been  made  by  Carter  [26]  who  reproves  result  (i)  with  the  changes 
indicated,  except  that  he  needs  to  assume  axial  symmetry  -  an  unfortunate 
restriction  which  no  one  has  yet  succeeded  in  removing. 

In  the  case  where  angular  momentum  is  nonvanishing,  established 
results  have  been  found  only  recently.  Carter  [27]  was  the  first  to 
provide  valuable  evidence  in  support  of  our  conjecture.  He  again  considers 
the  charge-free  case,  and  shows  that  all  asymptotically  flat,  stationary 
vacuum  fields,  satisfying  several  additional  (and  unwanted)  requirements 
including  axial  symmetry,  comprise  a  2-parameter  family,  where  one  of  the 
parameters  is  the  magnitude  of  the  angular  momentum.  Moreover,  only  the 
family  of  Kerr  solutions  contains  a  member  with  J  =  0  ,  i.e.  only  the 
Kerr  family  contains  the  Schwarzschild  subfamily. 

Very  recently  Hawking  [28]  has  proved  the  following  important 
result:  every  stationary  field  containing  an  event  horizon,  which  consists 

of  past  and  future  sheets  intersecting  in  a  2-space  with  the  topology  of  a 
2-sphere,  and  which  is  vacuum  in  some  neighbourhood  of  the  horizon,  is  eithe 
static  or  axisymmetric .  This  theorem  can  now  be  combined  with  the  earlier 
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theorems  of  Israel  (static)  or  Carter  (axisymmet ric)  to  "almost  prove" 

(we  still  have  several  unwanted  requirements)  that  the  field  must  be 
Schwarzschild  or  Kerr,  respectively. 

If  this  "combined  theorem"  eliminating  the  possibility  of  other 
2-parameter  families  distinct  from  the  Kerr  family,  can  be  reproved,  without 
the  unwanted  conditions,  it  will  establish  that  a  stationary  vacuum  black 
hole  is  a  Kerr  black  hole,  and  from  there,  it  seems  very  reasonable  that  the 
external  field  of  a  stationary  electrovac  black  hole  will  be  of  the  "charged 
Kerr"  type.  However,  at  present  these  remain  unverified  conjectures. 

§1.4  Summary  of  the  Thesis. 

The  previous  three  sections  have  provided  us  with  a  brief  survey 
of  the  overall  field  of  relativistic  gravitational  collapse.  The  purpose 
of  this  dissertation  is  to  present  several  separate  results  -  each  of 
which  deals  with  a  particular  phase  of  this  complex  problem. 

In  Chapter  II  we  present  some  fairly  elementary  results,  related 
to  the  instability  stage  of  spherical  collapse.  For  simplicity  we  deal 
with  spherical  shells.  A  solution  of  the  Einstein  equations  is  derived, 
representing  a  thin  spherical  shell  of  charged  fluid,  falling  in  a 
spherically  symmetric  field  due  to  mass  and  charge  at  its  centre.  No 
restrictions  are  placed  on  the  equation  of  state.  We  integrate  the  equations 
of  motion  to  find  the  law  of  conservation  of  total  energy,  and  we  use  it 
to  study  the  equilibrium  states  of  the  system,  and  their  stability 
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against,  collapse.  We  find,  under  reasonable  assumptions,  that,  given 
the  entropy  and  the  equation  of  state,  there  is  a  maximum  equilibrium 
mass,  and,  corresponding  to  it,  a  critical  radius,  inside  of  which 
instability  sets  in.  For  uncharged  bodies,  these  results  completely 
parallel,  and  serve  as  a  simple  illustration  of,  the  much  more  complica¬ 
ted  analyses  needed  for  fluid  spheres. 

Chapter  III  contains  studies  concerned  with  the  implosion  and 
horizon  stages  of  nearly-spherical  collapse.  Two  idealized  collapse 
models  (again  thin  shells) ,  involving  a  scalar  monopole  and  a  magnetic 
dipole,  are  considered,  treating  departures  from  sphericity  as  small 
perturbations.  Radiative  leakage  (largely  downwards  through  the  Schwarz- 
schild  horizon)  causes  the  externally  observable  asymmetries  to  decay  to 
zero  in  an  oscillatory  fashion.  These  results  have  significant  conse¬ 
quences  for  astrophysics;  they  imply  in  particular  that  a  black  hole 
cannot  be  a  source  of  synchrotron  radiation. 

In  Chapter  IV  we  deal  mainly  with  the  horizon  stage  of  highly 
nonspherical  collapse,  in  the  case  of  static  scalar  fields.  We  prove 
the  following  theorem.  Every  zero-mass  scalar  field  which  is  gravitation 
ally  coupled,  static  and  asymptotically  flat,  becomes  singular  at  a  simply 
connected  event  horizon.  In  the  special  case  where  the  gravitational 
coupling  of  the  scalar  energy  density  is  neglected,  the  solutions  are 
computed  explicitely.  Some  properties  of  the  singular  event  horizons  are 


also  discussed. 
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CHAPTER  II 

Gravitational  Instability  and  Collapse 

of  Charged  Fluid  Shells 


§2.1  Introduction. 

In  this  chapter  we  consider  the  motion  of  a  spherical  shell 
of  ideal  fluid,  with  charge  e  and  mass  m  ,  falling  in  the  external 
field  due  to  a  spherical  distribution  of  gravitational  mass  m^  and 
charge  e^  at  its  centre.  In  Newtonian  theory,  the  dynamics  of  this 
system  are  described  by 

T  +  U  +  (e2-e;L)2/2R  +  (e2-e1)e1/R  -  (n^-m^^R 

-  =  const»  » 

expressing  conservation  of  total  energy,  along  with  the  adiabatic 
condition 

dU  =  -P  d(4nR2)  , 

where  T  =  kinetic  energy  ,  U  =  internal  energy,  P  =  surface  pressure, 
e2  =  e+e^  ,  m2  =  m+m^  . 

It  is  our  aim  to  establish  the  relativistic  analogue  of  this, 


namely 


■ 


: 

12 


i  +  r2  =  a  +  b/r+  c/r2  , 

A  =  (m2-m1)2/M2  ,  B  =  mj-hn2  -  (n^-m.^)  (e^-e.^2) /M2  \  (2.1) 

2  2  2?  2  2  2 
C  =■  (e2  -ei  YlbYT  +  VT  lb  -  (e^+e^  )/2  , 

2 

where  A  ,  B  ,  C  are  functions  of  M  =  Mq  +  U  ,  and  dM  =  -Pd(47rR  )  . 

In  what  follows,  we  will  associate  with  the:  shell  three  different 
"masses":  Mq  (sum  of  the  rest  masses  of  the  constituent  particles 
of  the  shell),  the  total  proper  mass  M  (the  "bare  or  baryon  mass" 

Mq  plus  the  internal  thermal  energy  of  the  shell) ,  and  the  gravita¬ 
tional  mass  m  (the  total  energy  of  the  shell).  Mq  and  m  are  then 
constants  of  the  motion. 

Although  the  problem  of  gravitational  collapse  of  a  thin 
shell  is  a  hypothetical  one,  from  a  physical  point  of  view,  it  is  never¬ 
theless  useful,  because  in  its  crudest  overall  aspects,  the  dynamics  of 
stellar  collapse  turns  out  to  be  not  too  sensitive  to  the  detailed 
distribution  of  matter  within  the  star.  The  advantage  of  such  a  shell 
model  is  that  a  complete  solution  can  be  easily  found  in  an  explicit 
form  [29] . 


Various  special  cases  of  the  general  result  (2.1)  have  been 
obtained  previously.  Israel  [30]  has  derived  the  equation  of  motion  of 
a  spherical  shell  of  incoherent  dust  (m^=e^=e2=0  ;  zero  pressure)  and 
studied  its  collapse  [31] .  Using  a  similar  approach,  de  la  Cruz  and 
Israel  [32]  have  dealt  with  a  charged  shell  of  dust  in  the  field  due  to 
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a  concentration  of  mass  and  charge  near  its  centre  (zero  pressure) . 
Papapetrou  and  Hamoui  [33]  have  studied  uncharged  shells  of  ideal 
fluid  (m^=e^=e2=0)  and  their  motions,  assuming  a  particular  equation 
of  state  [34] .  Kuchar  [35]  has  considered  charged  fluid  shells 
(ml"ei=0)  ,  but  restricts  the  equation  of  state  to  be  polytropic.  We 
propose  to  obtain  (2.1)  without  placing  any  restriction  on  the  equation 
of  state. 


After  some  preliminaries  (§2.2),  the  equations  of  motion  for 
the  system  are  derived  (§2.3)  and  integrated  (§2.4),  giving  (2.1).  In 
§2.5,  we  deal  with  the  stability  of  an  uncharged  fluid  shell  (m^=e^= 
=e2=0)  in  equilibrium.  A  brief  treatment  of  stability  for  the  more 
complicated  original  system  is  given  in  §2.6.  At  several  points  we 
refer  to  Appendix  A,  which  contains  some  results  of  two-dimensional 
relativistic  gas  theory,  necessary  in  this  work. 


§2.2  Geometrical  and  Dynamical  Preliminaries. 

A  brief  sketch  of  the  geometry  and  dynamics  of  a  thin  shell  of 
ideal  fluid  is  offered  here  [36] . 


Let  £  be  a  timelike  hypersurface,  dividing  space-time  into 
two  4--dimensional  domains,  V+  and  V  ,  and  belonging  to  the  boundary 
of  both  [37] .  Suppose  are  intrinsic  coordinates  for  £  ,  with 


the  associated  tangent  base  vectors,  and  n  a  unit  spacelike 
normal  to  £  (directed  from  V_  to  V+)  .  For  an  arbitrary  vector  field 
A  ,  the  intrinsic  covariant  derivative  with  respect  to  is  defined  by 


' 

.a  '  a  .n  ff'Usvi  bn  .  r  ,)  £,S2 


■ 

" 
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Ai;  j  =  S(i)  •  -  Akrk> 


ij  * 


(2.2) 


where  the  Christoffel  symbols  r,  are  given  by 

k ,  i  j 


rk,ij  =  ~(k)  *  8~(i)/8^  * 


(2.3) 


Non-intrinsic  properties  enter  by  means  of  the  extrinsic  curvature 
3-tensor  of  £  ,  which  measures  the  variations  3n/9£1  of  the  unit 
normal.  Each  of  these  three  vectors  is  perpendicular  to  n  (hence 
tangent  to  [  ) ,  so  we  may  write 


j  • 

3n/3£  =  K.J  e 


i  ~(j) 


(2.4) 


which  leads  to 


K  —  e 

ij  ~(  j ) 


•  “  -5  •  8£(j)/351  =  -£  •  3£(1)/3CJ  “  K 


ji 


(2.5) 


In  general,  when  measured  with  respect  to  V+  and  V  ,  the  correspond¬ 
ing  components  of  this  tensor,  and  ,  may  or  may  not  be  equal. 

However,  in  order  that  £  be  (the  history  of)  a  surface  layer  [38],  we 
require  ^  K.^  .  It  then  follows,  provided 


T±j  =  K±j  -  K  -  [Ky  ] 


(2.6) 


is  nonvanishing,  that  J  is  the  history  of  a  thin  shell. 


From  (2.3)  and  (2.5)  we  obtain  the  Gauss-Weingarten  equations 


.  Ifiimron 
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n  +  r 


k 

ij  ~(k) 


(2.7) 


which,  with  (2.2)  leads  to 


3A/HJ  =  A1  ,  e , . . 
~  ;j  ~(i) 


-  A1  K . .  n 
ij  ~ 


(2.8) 


Taking  3/9£^  of  (2.7),  using  (2.4)  and  (2.7),  and  applying 
the  Ricci  commutation  relations  [39]  we  obtain  the  well-known  Gauss- 
Codazzi  equations 


Ra8y6  6(a)e(b)e(c)e(d)  =  Rabcd  "  Kac\d  +  Sc^d 


(2.9) 


R 


,a  ^ 


a8y6 


n  e(b)e(c)e(d)  =  \c;d  ^d; 


(2.10) 


Operating  on  (2.9)  and  (2.10)  with  g^Cga<^  and  g^ 
respectively,  and  using  the  fact  that 


be  3  y  By  By 
8  eCb)e(c)  =  8  -  n  “  1 


(2.11) 


we  find 


K  lK 
ab 


ab 


2  G  Q  n°nP 
aB 


(2.12) 


K  ,  -  K 
a  ;b  ;a 


r  a  nB 
G  0  6,  vll 
aB  (a) 


(2.13) 


3  r  ab 

where  R  is  the  intrinsic  3-curvature  invariant  of  l  .  K  =  g  Kab  , 

and  G  =  R  g  R  is  the  Einstein  tensor. 

a£  aB  2  °aB 
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Let  u  be  a  unit  vector,  tangent  to  £  ,  representing  the 
velocity  of  matter  in  the  shell.  Then,  from  (2.8),  the  4-acceleration 
of  an  element  of  the  shell,  as  measured  in  V+  and  V_  ,  is  [40] 


Su°76t 


+ 

—  a 


=  e 


u 


(i)  ;j 


•  j  • 

j  aT,  i  i 
,uJ  -  n  K^u  uJ 


(2.14) 


Since  any  motions  of  the  shell  are  in  a  direction  normal  to  the  shell, 
we  restrict  our  interest  to  the  normal  components  of  (2.14),  namely 


+  .  .  + 

n  6ua/6x  =  -u1u'*K.  . 
a  ij 


from  which,  by  adding  and  subtracting,  we  obtain 


n  6ua/6x 
a 


+ 


+  n  6ua/Sx 
a 


i  1  ~ 

-  2u  uJKij  , 


(2.15) 


n  5ua/6x 
a 


-  n  6ua/6x 
a 


=  -  u^"u^  (from  (2.6)).  (2.16) 


Turning  to  dynamical  considerations,  the  components  of  the 
surface  energy  3-tensor  are  defined  by  the  "Lanczos  equations", 


Yij  -  8ijY  -  -8TlSij  • 


(2.17) 


or  equivalently 


-8tt(S 


ij 


1 

2  8ij 


S) 


(2.18) 


We  also  have  Einstein's  field  equations. 
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(2.19) 


satisfied  in  the  regions  V_  and  V+  ,  exterior  to  the  shell.  Using 
(2.17)  and  (2.19),  the  jumps  of  (2,12)  and  (2.13)  across  £  can  be 
put  in  the  form 


SiJK  =  [T  n°n8]  ,  (2.20) 

ij  aft 

Si;  j  ‘  -[TaSe(i)nSl  •  (2’21> 

For  a  shell  of  ideal  fluid,  the  intrinsic  surface  energy  tensor  has 
the  form 


Sij  =  (a+P)  u^u^  +  Pg1^  ,  =  ~  1  >  (2.22) 

where  P  is  the  surface  pressure,  and  a  is  the  surface  energy  density 
of  the  fluid.  It  is  easily  shown  from  (2.20)  and  (2.22)  that 

-2u1u^K  =  2 (a+P)-1  (PK-[T  nan6]}  , 

J- j  CXp 

and  from  (2.18)  and  (2.22)  that 


-u^u^y^j  =  87r(P+a/2) 


Hence  (2,15)  and  (2.16)  can  be  written 


n  6ua/6i 
a 


+  n  6ua/6x 
a 


2 (a+P)  1  {PK-[Tagnane]>  , 


(2.23) 


^  ■  ~  bn  r .  )  i  bna 
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n  6ua/6x 
a 


+ 


-  n  6ua/6x 
a 


=  8tt(P+o/2)  . 


(2.24) 


Note  that  these  results  agree  with  those  of  de  la  Cruz  and  Israel  [41] 
for  zero  pressure,  and  are  identical  with  Kuchar's  results  [42], 


For  future  reference,  it  is  easily  derived  from  (2,21)  and 
(2.22)  that 


(cru^)  .  .  +  Pu^ 

» J  » J 


'TaBe(i/]  U±  “  tT<*uVl 


(2.25) 


§2.3  Charged  Spherical  Shell  of  Ideal  Fluid  in  a  Spherisymmetric 
Electrovac  Field. 


Consider  a  charged  shell  of  ideal  fluid,  falling  in  the  field 
due  to  a  spherically  symmetric  distribution  of  mass  m^  and  charge 
e^  ,  near  its  centre.  Let  r  =  R(t)  be  the  equation  of  £  .  Then  its 
intrinsic  metric  is  given  by 


(ds2)^  =  (R(t)}2  dft2  -  dx2 


(dft2  =  d02  +  sin20d(f>2) 


(2.26) 


where  £  =  (0,<]>,x)  and  x  is  proper  time  along  the  streamlines 
0,4>  =  constant.  In  both  V+  and  V  ,  the  line  element  is  reducible  to 
the  Reissner-Nordstrom  metric,  by  extension  of  Birkhoff’s  theorem  [43], 
Hence  we  have 


(ds2)_  =  {f _(r) }  1dr2  +  r2dft2  -  f_(r)dt2 
(ds2)+  =  {f+(r)}_1dr2  +  r2dfi2  -  f+(r)dt2 


(r  <  R(x) ) 


(2.27) 


(r  >  R(x) ) 


(2.28) 


3Bri3  (££.£) 
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where 


f _(r)  =  l-2m^/r  +  e2/r2  ,  f+(r)  =  l-2m2/r  +  e2/r2  . 


(2.29) 


It  is  easily  seen,  by  proper  choice  of  t  (x)  ,  t+(x)  on  £  ,  that 
the  metrics  induced  by  (2.27)  and  (2,28)  on  £  ,  agree  with  (2.26). 

It  can  be  verified,  by  intrinsically  differentiating 


u°u  =  -1  ,  with  R  =  dR/dx  and  t.  =  dt,/dx  ,  that  [44] 
a  +  + 


0  ■  u  6ua/6x 
a 


=  f  (R) _1  R  6  R/<5x  -  {f , (R)  +  R2}2  6  t  / 6x  .  (2.30) 

+  +  + 


We  now  use  (2.30)  to  eliminate  6t_)_/6x  from 


n  6ua/6x 
a 


JL 

2,2 


=  (f  (R)  +  R^}Z  f  (R)  1  6  R/6x  -  R6  t+/6x 


and  since 


6R/6x  -  R  +  j  df+(R)/dR  , 


we  finally  obtain,  from  (2.29), 


n  Su  /Sx 
a 


1 

a/"  =  {f+(R)  +  R2}  2  (R  +  m2/R2  -  e2/R3}  . 


(2.31) 


ot 

The  analogous  expression  for  n  6u  / 6x 


a 


is 


n  5ua/6x 
a 


-  (f  (R)  +  R  } 


1 

2 


(R  +  m^/R^ 


e|/R3} 


(2.32) 


to  soiorto  laqo^q  .nase  x  ni  3 1 


s  i  Bi'  £)  (  ;,£)  <>2W  W-  »W 


-C 


20  - 


The  energy  tensor  associated  with  the  Reissner-Nordstrom 
metric  (2.27)  is  given  by 

^  1  2  3  2  A 

-T^  =  -T^  =  =  e^/8Trr  (other  components  zero),  (2.33) 


with  e^  replaced  by  e ^  for  the  metric  (2.28).  From  (2.33)  it  is 
easily  shown  that 


T^u°n 

a  3 


=  0 


(2.34) 


[T^nane]  =  -(e2  -  e2)/87rR4 


(2.35) 


Using  (2,25)  and  (2.34)  it  follows  that 


(auJ)  =  -Pu^  .  (2.36) 

» J  J  J 

The  left  side  of  (2.36)  is  just  the  rate  of  increase  of  surface  energy, 
while  the  right  side  is  equal  to  minus  the  rate  of  work  done  by  the 
pressure  in  expanding  the  shell.  Hence,  unlike  the  de  la  Cruz  and 
Israel  result  for  dust  [45],  the  proper  mass  of  an  element  of  fluid 
under  surface  pressure  P  changes  with  time. 


We  can  now  write  (2,23)  and  (2,24),  using  (2.31),  (2.32)  and 
(2.35),  in  the  form 


,{f+(R)  +  R  } 


1 

2 


{R+m2/R2  -  e2/R3}  +  (f_(R)  +  R2 } 


1 

2 


*'  2  2  3 

(R+n^/R  -e1/RJ} 


=  2(a+P)"1(PK+  (e2-e2)/87iR4}  , 


(2.37) 


'  i  n-JI3 
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(f+(R)  +  R  } 


1 

2 


{R+m2/R2  -  e2/R3}  -  {f  (R)  +  R2} 


1 

2 


,2  2  3 

{R+m^/R  -  e^/R  } 


=  87r(P+a/2) 


(2.38) 


We  call  (2.37)  and  (2.38)  the  equations  of  motion  for  the  system.  In 

•2  1/2 

order  to  simplify  the  notation,  we  define  F+  =  (f+(R)  +  R  }  and 

F_  =  (f_(R)  +  p2 } 1 / 2  ,  Then  (2.37)  and  (2.38)  become  simply 

{F+  +  F_ }/R  =  2(q+P)“1  (PK  +  (e2-e2)  8ttR4}  ,  (2.39) 

{F+  -  F_}/R  =  8tt (P+a/2)  .  (2.40) 


§2.4  Integration  of  the  Equations  of  Motion. 


A  first  integral  of  (2.40)  is  obtained  by  noting  that,  for 
the  intrinsic  coordinates  of  (2.26),  we  have  u1  =  (0,0,1)  ,  therefore 
u^  =  2R/R  .  From  (2.36),  it  then  follows  that 


(P+a/2)R  =  -  -j  d(aR)/dx 


(2.41) 


•  • 

It  is  worthy  pointing  out  that  (2.41)  is  just  M  =  -PA  ,  where  A  = 

2 

area  =  4ttR  ,  and  M  is  the  total  proper  mass  of  the  shell,  defined  by 


M 


/  r.2 

4ttR  a 


(2.42) 


Using  (2.41),  (2.40)  can  be  written  equivalently  as 


. 
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F+  -  F_  =  -4mi(aR)  /d  x  ,  (2.43) 

which  can  be  integrated  to  give  (from  (2.42)) 

F+  -  F_  =  -  M/R  -  C  ,  (2.44) 


where  C  is  a  constant  of  integration.  Multiplying  (2.44)  by 
(F  +  F  )  and  using  (2.42)  leads  to 

F+  +  F_  =  2(m2-m1)/R(47TR0+C)  -  (e2-e2)  /R2(47tRct+C)  .  (2.45) 

From  (2.5),  a  long  but  straightforward  calculation  gives 

K  -  Kijglj  =  {F+  +  F_ }/R  +  j  {F+  +  F_}/R  .  (2.46) 

Putting  (2.46)  into  (2.39)  and  simplifying,  yields 

9  9  / 

{F+  +  F_}  =  (2PR/aR)  (F+  +  F_)  +  (e2  -  e  )R/4ttR  a  .  (2.47) 

Substituting  (2.45)  and  its  time  derivative  into  (2.47),  and  comparing 
both  sides,  with  the  aid  of  (2.41),  it  follows  that  C  must  vanish  in 
order  that  (2.45)  be  compatible  with  (2.39).  Hence,  using  (2.42), 
equations  (2.44)  and  (2.45)  become 

1  1 

M  =  F_R  -  F+R  =  {R2(l+R2)  -  2mxR  +  e2}  2  -  (R2(l+R2)  -  2m2R  +  e2}2  , 

(2.48) 


2  2 , 


F_R  +  F+R  =  {2R(m2-m1)  -  (e2-e  )}/M 


(2.49) 
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Differentiation  of  (2.49),  with  the  help  of  (2.41),  will  give  (2.47), 
showing  that,  as  expected,  (2.49)  is  an  integral  of  (2.39).  Equations 
(2.48)  and  (2,49)  can  be  combined,  by  adding  and  squaring,  to  give 
the  final  form,  (2.1). 

At  this  point,  several  remarks  should  be  made  concerning  the 
integral  of  motion. 

(i)  This  integral  (2.1)  along  with  (2.41)  (i.e.  along  with 
2 

dM  =  -Pd(4frR  ))  sums  up  the  information  contained  in  all  the  previous 
equations  related  to  the  motion  of  the  system. 

(ii)  Our  result  is  precisely  that  obtained  by  de  la  Cruz  and 
Israel  [46]  for  zero  pressure,  except  that,  in  our  case,  the  total 
proper  mass  M  is  not  constant. 

(iii)  In  the  special  case  m^  =  e^  =  0  ,  considered  by  Kuchar  [47], 

(2.1)  can  be  written  in  the  equivalent  form 

1 

m  =  M(l+R2) 2  -  (M2-e2)/2R  ,  (2.50) 

expressing  the  conservation  of  total  energy  m  of  the  shell.  Our 
integral,  then,  is  a  generalization  of  this  conservation  equation,  to 
non-zero  m^  and  e^  . 

(iv)  The  integral  (2.1)  is  also  a  generalization  of  Kuchar* s 
result  in  a  different  sense.  In  obtaining  (2.50),  Kuchar  makes  the 
assumption  that  the  equation  of  state  is  polytropic,  i.e.  that  the 
adiabatic  exponent,  defined  by 


.sense  3/rsisinb  s  nl  jLub si 


Y  =  -  d£nP/d£nA 


» 


(2.51) 


is  constant.  In  our  considerations,  however ,  no  mention  has  been  made 
of  any  restrictions  on  the  equation  of  state,  and,  in  fact,  none  need 
be  made . 


In  Section  2.6  we  consider  the  stability  of  this  system,  and 
we  make  use  of  definition  (2.51),  keeping  in  mind  that  y  is  not 
constant,  but  is  itself  a  function  of  the  entropy  and  the  density  [48] . 


§2.5  Stability  of  an  Uncharged  Fluid  Shell. 

In  this  section,  we  avoid  several  involved  computations, and 
consider  the  reasonably  simple  case  of  an  uncharged  fluid  shell 
(m^  =  e^  =  =  0)  .  We  want  to  study  the  stability  of  such  a  system 

against  collapse. 

The  integral  (2.1)  can  be  simplified,  in  this  case,  to 
1 

9  9  9 

m  =  M( 1+R  )  -  M  /2R  =  total  energy  ,  (2.52) 

which,  along  with 


dM  =  -  Pd(47rR2)  ,  (2.53) 

and  the  equation  of  state,  completely  describes  the  motion  of  the 
shell.  Consider  a  momentarily  static  configuration,  that  is,  a  shell 
for  which,  at  some  instant,  R  =  0  .  Taking  d/d x  of  (2,52),  and 
using  (2.53),  it  is  immediate  that 


i  ft 
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(MR)  ^  =  8ttR(1-M/R)  (P-r)  ,  (2.54) 

R=0 

where  r  is  defined  by 

T  =  M2/16ttR2(R-M)  ,  (2.55) 

representing  the  gravitational  self -attraction  of  the  shell  [49].  If 

the  equations  of  state  P  *  P(oq)  ,  a  =  a(ao)  are  given  (cq  is  the 

"bare  mass  density"),  then  the  configuration  and  initial  acceleration 

of  any  momentarily  static  shell  are  completely  determined  by  any  two 

pieces  of  data.  For  example,  suppose  we  are  given  the  number  of 

particles  (i.e.  given  M  )  and  the  radius  R  .  Then  a  follows 

2 

from  oq  =  Mo/4ttR  ,  and  hence  P  =  P(oq)  ,  a  =  a(aQ)  •  M  ,  m  , 

and  r  are  found  from  (2.42),  (2.52)  and  (2.55)  respectively.  Finally 

R  follows  from  (2.54).  Hence,  to  any  momentarily  static  configuration 

of  a  given  shell  (of  fixed  M  ) ,  we  can  associate  a  potential  function 

m(R,M  )  : 
o 

m#  =  M  -  M2/2R  =  M  +  U  -  M2/2R  .  (2.56) 

R=0,M  fixed 
o 

In  a  quasistatic  displacement,  the  change  in  this  potential  function 
can  be  written 

dm.  -  -  8ttR(1-M/R)  (P-r)dR  , 

R=0,M  fixed 
o 

1  /2 

and  since  1-M/R  =  (l-2m/R)  from  (2.56),  we  have 
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1 

dm.  =  -(l-2m/R) 2  (P-r)d(47TR2)  .  (2.57) 

R=0,M  fixed 
o 

This  can  be  thought  of  as  the  energy  change,  measured  by  an  observer 

2 

at  infinity.  The  term  (P-r)d(47rR  )  is  the  local  measure  of  work 
done,  AE  ,  and  it  is  decreased  by  ~AEm/R  ,  the  amount  of  energy  lost 
in  pushing  AE  up  to  the  distant  observer. 

Suppose  now  that  the  shell  is  in  equilibrium.  Then 

R  »  R  =  0  ,  (2.58) 


hence  from  (2.54) 


p  =  r 


From  (2.57),  this  is  seen  to  be  equivalent  to 


(2.59) 


(dm/dR).  =  0  (2.60) 

R=0,M  fixed 
o 

% 

i.e.  the  potential  energy  is  stationary  for  equilibrium. 

As  a  result  of  this  additional  condition,  we  have  that, 

given  the  functions  P(a  )  »  tf(aQ)  »  an  equilibrium  state  is  completely 

determined  by  a  single  datum.  Under  the  assumption  of  equilibrium, 

2 

P  =  T  ,  equation  (2.55),  with  M  =  4nR  a  ,  gives  R  as  a  function  of 
P  ,  a  : 


R  =  P/iTa(a+  4P) 


(2.61) 


' 


•V 


. 
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Similarly 


M  =  4P^/7ra(cH-  4P)^ 


(2.62) 


and 

m  =  4P^(a+2P) /Tra(a+4P)  ^  (2.63) 

follow  from  (2.42)  and  (2.56). 

Consider  the  class  of  all  equilibrium  states,  each  of  which 

is  specified  by  a  particular  equation  of  state,  and  one  piece  of  data 

(say  Mq)  .  In  the  non- relativistic  limit  (P  <<  a  )  ,  the  equation 

2 

of  state  is  given  by  P  a  a  (see  Appendix  A),  hence,  from  (2.63), 

lim  m(P,a(P))  =  0  .  Similarly  for  the  ultra- relativistic  limit 

P-*0 

(P  >>  a  )  ,  we  have  P  a  a  (see  Appendix  A)  and  lim  m(P,a(P))  =  0  . 

°  p-*x> 

We  conclude  that  there  must  be  at  least  one  maximum  for  m(P,a)  .  If 
we  suppose,  for  simplicity,  that  there  is  exactly  one  maximum,  we  can 
represent  the  class  of  all  equilibrium  states  P  =  r  ,  by  the  solid 
curve  in  Figure  1.  For  a  given  m  ,  then,  there  may  be  two  or  more 
equilibrium  states.  Presumably,  not  all  of  these  will  be  stable 
against  collapse. 

Let  us  look  at  the  stability  of  a  particular  shell  of  fixed 
Mq  ,  whose  equilibrium  configuration  is  represented  by  the  point  E 
(Figure  1).  The  assignment  of  Mq  completely  determines  m^,  ,  R^,  , 
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The  (m,P)  plane.  Given  the  equation  of  state  P  =  P(a)  ,  each  point  in  this  plane 
represents  a  momentarily  static  configuration  (MSC)  of  some  shell.  All  MSC's  accessible 
to  a  shell  containing  a  definite  number  of  baryons  (fixed  M  )  lie  on  a  curve  in  the  (m,P) 
plane  (the  dashed  curve  corresponds  to  one  particular  value  of  MQ)  .  The  curve  P  =  p 
represents  all  possible  equilibrium  states. 
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Pg  .  We  now  displace  the  shell  quasist atically  to  A  ,  as  in  (2.57), 
keeping  Mq  fixed.  Since  work  is  done,  mA  ^  .  At  A  we 

release  the  shell  to  move  freely  (m  =  m^  =  constant) .  In  order  that 
the  system  be  stable,  it  is  necessary  and  sufficient  that  R  is 
opposite  in  sign  to  ,  thus  causing  the  shell  to  move  from 

A  to  B  ,  turn  around,  and  continue  to  oscillate  about  E'  .  Hence, 
from  (2.54),  (P-T)^  must  be  opposite  in  sign  to  R^  -  R^  .  Therefore, 

a  necessary  and  sufficient  condition  that  the  equilibrium  at  E  be 
stable  is  that 

(d/ dR) (P-r)  <0  ,  (R  =  R  )  .  (2.64) 

R=0,M  fixed 
o 

From  (2.57)  and  (2.60),  we  see  that  (2.64)  is  equivalent  to 

(d2m/dR2)#  >0  ,  (R  =  R  )  ,  (2.65) 

R=0,M  fixed 
o 

i.e.  the  equilibrium  at  E  is  stable  if  and  only  if  the  potential 
function  m(R^,Mo)  is  a  minimum.  In  the  remainder  of  this  section,  we 
will  understand  that  dm/dR  etc.  are  taken  with  R  =  0  ,  Mq  fixed, 
although  we  no  longer  write  these  conditions,  for  brevity. 

We  now  derive  an  expression  for  the  critical  radius,  R^,  ,  at 

which  stable  equilibrium  configurations  become  unstable.  Taking  d/dR 
2  2 

and  d  /dR  of  (2.56)  at  R  =  R^  ,  conditions  (2.60)  and  (2.65)  become 

dm/dR  =  (dM/dR)E(l-M/RE)  +  M2/2R2  =  0  ,  (2.66) 
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d2m/dR2  =  (d2M/dR2)E(l-M/RE)  -  (dM/dR)g((dM/ dR)£  -  2M/Rg)/Rg  - 

-  M2/Re  >  0  (2.67) 

From  (2.51)  and  (2.53),  it  is  not  difficult  to  show  that 

(d2M/dR2)E  =  (dM/dR)E(l-2y)/RE  .  (2.68) 

Substituting  (2.68)  into  (2.67)  and  simplifying,  by  means  of  (2.66), 
the  stability  condition  becomes 


y  -  3/2  >  (M/(Re-M)}  (M/4(Re-M)  +  1}  .  (2.69) 


Using  (2.56),  (2.69)  can  be  written 

1  1 

Y  -  |  >  {-1  +  Re/  (R^-2REm)2}  {i(iy  (R^^m)2)  +  |}  .  (2.70) 


To  obtain  an  approximation,  we  assume  Rg  >>  m  and  neglect  second 
order  terms  in  m/Rg  .  Then  (2.70)  becomes  simply  the  condition  that 
the  equilibrium  at  E  is  stable  if  and  only  if 


Rg  >  m/(y  -  -j)  ~  Rc  ,  <  y  _<  2 


(see  Appendix  A)  .  (2.71) 


It  is  interesting  to  note  the  similarity  of  this  result  and  that  of 
Chandrasekhar  [50]  for  fluid  spheres,  namely  that  stable  equilibrium 
occurs  for 


Rg  >  2  Km/  (y  _'f)~Rc  >  <  Y  £  >  K~1  . 


' 
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We  conclude  this  section  by  showing  that  equilibrium  states 
such  as  E  (where  P  <  P^)  are  stable  against  collapse,  while  config¬ 
urations  such  as  F  (P  >  P^)  are  unstable.  Since  dm/dP  =  (dm/dR)  * 
(dR/dP)  ,  we  see  from  (2.59)  that  the  =  constant  curve  (broken  line 
in  Figure  1)  has  extrema  at  E  and  F  .  It  also  follows  from  (2.60) 
that 

d2m/dP2  =  (d2m/dR2)(dR./dP)2  ,  (2.72) 

hence  from  (2.65)  the  equilibrium  at  E  is  stable  if  and  only  if  that 

curve  has  a  minimum  there.  Since  no  extremum  occurs  between  E  and  F  , 

it  is  immediate  that  stability  (instability)  at  E  implies  instability 

(stability)  at  F  .  This  argument  shows  that  the  critical  radius  R^ 

must  correspond  to  the  equilibrium  state  at  C  .  That  the  equilibrium 

is,  in  fact,  stable  at  E  ,  can  be  shown  by  considering  (2.71)  in  the 

non-relativistic  limit  (y  +  2)  .  In  this  limit,  as  E  approaches  the 

2 

origin  along  P  =  T  ,  we  have  P  a  a  ,  hence  lim  m(P,a(P))  =  0  , 

P+0 

while  lim  R(P,a(P))  is  finite,  from  (2.61).  Therefore  (2.71)  is 
P+0 

satisfied  near  the  origin  and  we  have  stability  there.  Consider  again 

a  quasistatic  displacement  from  E  (near  the  origin)  to  A  along 

M  =  constant.  If  we  make  the  reasonable  assumption  that  P  is  a  mono- 
o 

tone  function  of  a  ,  then  dM  =  0  implies 

o  o 

(dR/dP)  =  -  (R/2a  )  (da  /dP)  <  0  . 

o  o 


Thus,  the  increase  in  P  from  E  to  A  is  coincident  with  a  contraction 


' 
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of  the  shell.  This  fact,  together  with  the  stability  at  E  ,  is 
sufficient  to  show  P  >  r  throughout  the  region  under  the  P  =  r 
curve.  It  then  follows  that  the  equilibrium  states  are  stable  for 
P  <  Pc  and  unstable  for  P  >  Pc  . 


§2.6  Stability  of  the  Charged  System. 


We  return  now  to  the  original  system,  and  derive  a  condition 

for  stability  against  collapse.  We  consider  only  configurations  which 

are  outside  the  upper  Nordstrom  horizon,  i.e.  we  assume  [51] 

1  1 

•p  t  2  2,2  ^  f  2  2,2 

R  2.  m2  +  (m2  “  e2^  >  +  (m^  -  e^) 


It  is  convenient  to  deal  with  the  integrals  (2.48)  and  (2.49)  rather 
than  (2.1).  For  an  equilibrium  configuration  at  R  =  Rq  we  must  have 
(2.58)  satisfied  there.  Then  one  equilibrium  condition  can  be  written, 
from  (2.48)  and  (2.58),  as 


M=  (F_R-  F+R)r  =  {(R  >2  . 

o  o  1 


2  2?  2 
(m l-epr  -  {(Ro-m2)Z 


,  2  2,  2 

(m2-e2) } 


(2.73) 


Taking  d/dx  of  (2.48)  and  using  (2.58),  a  second  equilibrium  condition 
is  found  to  be 

(dM/dR)R  =  (Ro-m1)(F_R)'1  -  (R^)  (F+R) ^  . 

o  o  o 


(2.74) 
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Equivalently,  if  we  write  (2.49),  from  (2.48),  in  the  form 

m  =  n^-m^  =  F_M  -  {M^  -  (e^-e^) }/2R  =  constant  ,  (2.75) 

and  follow  an  argument  identical  to  Section  2.5,  we  see  that 

(P-r)R  =  0  ,  (2.76) 

o 

in  order  that  equilibrium  be  maintained.  In  this  case,  r  represents 
the  gravitational  self-attraction,  together  with  the  electrostatic 
forces,  and  is  given  by 

r  =  {M2  -  (e2-e2)  +  2M(m  R  -e2)(F  R)I1}/16ttR2  {(F  R)  -M}  . 

R  21  lol-R  o-R 

o  o  o 

(2.77) 

As  in  the  previous  section,  the  equilibrium  will  be  stable,  provided 

(d/dR)(P-r)>  <  0  ,  (R  =  R  )  (2.78) 

R=0,M  fixed 

If  we  now  substitute  for  P  from  (2.53),  and  for  r  from  (2.77), 
and  evaluate  the  left  side  of  (2.78),  the  resulting  complicated  expres¬ 
sion  can  be  reduced,  with  the  aid  of  (2.73)  and  (2.74),  to 

(d2M/dR2)  +  {(m2-e2)(F_R)‘3  -  (m2-e2) (F+R)~3 }  >  0  .  (2.79) 

o  o  o 

From  (2.68)  and  (2.74)  we  find 

(d2M/dR2)R  =  {(Ro-m1)(F_R)’1  -  (R^m.,)  (F+R)'L }  (1-2y) /Rq  .  (2.80) 

o  o  o 
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Substitution  of  (2.80)  into  (2.79)  leads  to  the  final  form 

{(Vml)(F-R)R1  "  (Vm2)(F+R)R1}  (1'2V)/R0  >  (m^XF  R)"3  - 

o  o  o 

-  (m^-e^)(F_R)^3  .  (2.81) 

o 

This  form  of  the  stability  condition  is  a  function  of  the  equilibrium 

3 

radius  Rq  ,  alone,  since  m^  ,  m2  ,  ,  e^  are  fixed.  Since  y  >  y  , 

2 

and  (dM/dR)R  is  negative  (from  dM  =  -Pd(4TrR  )  ,  P  >  0)  ,  it  is 
o 

clear,  from  (2.74),  that  the  left  side  of  (2.81)  is  positive.  Thus, 

2  2  1/2 

instability  for  some  Rq  >  m2  +  (m2~e2^  can  occur>  only  when  the 

right  side  of  (2.81)  is  positive.  We  now  show  that  this  is  always 
the  case. 


The  obvious  restriction  on  the  equation  of  state  that  P  >  0  , 
becomes,  for  equilibrium,  the  condition  that  r  >  0  .  From  (2.77), 
this  restriction  can  be  reduced  to  the  simple  condition 


(m2i-e2)/(Vmi)2  <  (m2-e2)  /  (RQ-m2) 2 


(2.82) 


2  2 

This  is  a  generalization  of  the  condition  m  -  e  >  0  ,  found  by 
Kuchar  [52].  Using  (2.82),  it  is  straightforward  to  show  that  the  right 
side  of  (2.81)  is  always  positive.  We  therefore  see,  under  the  physi¬ 
cally  reasonable  assumptions  |e^|  c  m^  ,  | e^ I  <  ,  that  instability 

can  occur  outside  the  upper  Nordstrom  sphere.  To  show  that,  in  fact. 


it  does  occur,  it  is  sufficient  to  prove  the  existence  of  a  critical 


2  2  1/2 

radius  Rq  >  m2  +  (m2~e2^ 


If  we  evaluate  both  sides  of  (2.81)  at 


2  2  1/2 

the  Nordstrom  sphere  R  =  +  (m^-e^)  ,  we  find  that  condition  (2.81) 


Ji/to  no:.3ib;  )o  art  ,  J.  i  ,r  05  s  acosvf 
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is  violated,  indicating  instability.  Similarly,  it  is  trivial  to  show 

that,  as  R  ■+  00  ,  the  condition  is  satisfied,  hence  instability  sets 

2  2  1/2 

in  somewhere  outside  R  =  +  (m2~e2) 

An  approach  which  is  perhaps  more  instructive,  is  to  assume 
Ro  >>  m 2  ,  and  find  an  approximating  stability  condition,  which  is 
more  informative  than  the  exact  relation.  Neglecting  second  order  terms 
in  m2/R0  >  we  can  write  (2.75),  with  R  =  0  ,  as 

m  =  M  +  \  { (m^-e^) / (R  -m  )  -  (m9-e9)/(R  -m9)  }  , 

A  II  Ol  A  A  O  A 

and  hence  an  approximating  condition  can  be  found,  from  (2.65),  to  be 


R 

o 


> 


2Y-1 

2y-3 


2  2 

m1(m1-e1) 


2  2 
(m^-ep 


m2(m2"e2 > 


2  2 

(m^ep 


(2.83) 


for  stability.  Under  the  assumption  Rq  >>  >  (2.82)  becomes  simply 

2  2  2  2 

m^-e^  <  m2~e2  »  r.e.  the  denominator  of  (2.83)  is  positive.  Hence 


R 

o 


2X1 

2Y-3 


» 


for  stability.  In  the  case  of  particular  interest,  an  effectively  "hot" 

3 

gas  exerting  pressure,  we  have  y  +  and  therefore,  only  very  diffuse 
equilibrium  states  are  stable. 
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CHAPTER  III 

Gravitational  Collapse  with  Asymmetries 


§3.1  Introduction 


Every  static  nonspherical  perturbation  of  Schwarzschild ' s 
exterior  field,  due  to  gravitational  or  electromagnetic  sources  within 
the  stationary  lightlike  surface  gQQ  =  0  [53]  becomes  singular  on  this 

surface  [54].  The  same  is  true  for  all  static  zero-mass  scalar  pertur¬ 
bations  [55]  -  spherical  as  well  as  nonspherical.  Further,  it  appears 
from  the  discussion  in  §1.3,  that  stationary  perturbations  of  Kerr's 
rotating  solution  have  a  similar  property  [56].  Assuming  these  results 
to  be  applicable  to  the  asymptotically  stationary  exterior  field  of  a 
collapsing  star,  one  is  led  to  the  conjecture  that  all  externally  detect¬ 
able  asymmetries,  including  scalar  and  magnetic  fields,  must  somehow  decay, 
leaving  behind  Schwarzschild ' s  vacuum  field  (or,  in  the  case  of  nonvanish¬ 
ing  charge  and  angular  momentum,  the  "charged  Kerr"  field)  as  the  sole 
external  manifestation  of  the  collapsed  object. 

To  examine  these  questions,  we  carry  out  a  dynamical  analysis  of 
two  idealized  collapse  models,  one  involving  a  scalar  monopole,  the  other 
a  magnetic  dipole.  Our  results  support  the  foregoing  conjecture  and  reveal 
the  decay  mechanism  to  be  a  rapid  radiative  leakage  of  the  perturbing  field, 
largely  downwards  through  the  event  horizon. 


We  cast  the  Schwarzschild  metric  into  the  form 
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2  2  2 
ds  =  a  dxdy  +  r  d^  , 


a  =  1  - 


(3.1) 


where  the  retarded  and  advanced  time  coordinates  -x,y  are  related  to 
the  standard  Schwarzschild  coordinates  by 


\ 


x  =  (r-1)  +  £n(r-l)  -  t  , 
y  =  (r-1)  +  £n(r-l)  +  t 


(3.2) 


Lengths  are  measured  in  units  of  the  Schwarzschild  radius  2m  =  1  . 


Both  of  our  models  can  be  considered  as  linearly  perturbed 
variations  of  the  following  basic  situation  (Figure  2).  A  thin  hollow 
spherical  shell  of  mass  m  =  y  is  initially  static  with  radius 
r  =  Rq  >>  1  [57].  At  time  t  =  —  |  x^  =  -(R-1)  -  £n(RQ-l),  it  suddenly 
begins  to  collapse  with  the  speed  of  light  [58].  This  model,  adopted  for 
mathematical  simplicity,  is  highly  artificial  from  an  astrophysical  point 
of  view,  but  does  not  violate  any  of  the  principles  of  relativity  theory. 
Moreover,  our  main  interest  is  in  the  asymptotic  behaviour  of  the  exterior 
field  as  t  -*■  00  ,  and  we  do  not  expect  this  to  depend  too  sensitively  on 
the  precise  structure  of  the  source  or  the  initial  conditions. 

The  perturbations  we  intend  to  impose  on  this  idealized  collapse 
picture  are  of  two  types  -  one  a  static  spherically  symmetric  scalar  mono¬ 
pole,  the  other  a  static  axi-symmetric  magnetic  dipole.  Both,  however,  are 
located  at  the  centre  of  the  collapsing  shell,  and  both  are  assumed  to  be 


weak. 
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Space-time  diagram  for  collapsing  shell  model. 
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Since  news  of  the  onset  of  collapse  cannot  reach  the  interior 
ahead  of  the  shell  itself,  the  initial  static  interior  field  (Region  I 
in  Figure  2)  remains  unchanged  in  both  models.  The  exterior  field,  how¬ 
ever,  becomes  time-dependent  after  passage  through  a  shock  front  at 
x  =  xq  .  Our  problem  is  thus  to  find  the  perturbing  field  in  the  time- 
dependent  region  x  _<  xq  ,  y  ^  0  (Region  III  in  Figure  2),  given  Cauchy 

initial  data  on  the  characteristics  y=0,x=x  >>1.  This  we  will 

o 

do  by  numerical  integration.  (It  is  unnecessary  to  distinguish,  at  this 
level  of  approximation,  between  null  hypersurfaces  of  the  Schwarzschild 
background  field  and  of  the  perturbed  field). 

In  Section  3.2  we  formulate  the  characteristic  initial  value 
problem  for  the  collapse  model  just  described,  in  the  case  of  the  scalar 
monopole  perturbation.  Section  3.3  contains  the  corresponding  derivation 
for  the  magnetic  dipole  case.  It  is  seen  that  the  two  problems  can  be  put 
in  very  similar  form  -  a  convenient  fact  when  it  comes  to  integrating  numer¬ 
ically.  In  Section  3.4  the  results  of  the  integrations  for  these  two  models 
are  presented,  including  a  discussion  of  where  the  external  scalar  or 
magnetic  field  energy  goes.  Finally,  Section  3.5  gives  qualitative  inter¬ 
pretation  and  some  concluding  remarks. 

§3.2  Scalar  Monopole  Model 

For  our  first  model  we  consider  a  static  spherically  symmetric 
zero-mass  scalar  monopole,  located  at  the  centre  of  the  shell.  It  is 
assumed  that  this  source  is  weak  enough  that  any  gravitational  effects  of 
the  scalar  energy  density  can  be  neglected  for  r  >_  1  (i.e.  the  scalar 
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field  is  not  gravitationally  coupled). 


In  order  to  find  the  perturbing  scalar  field  in  the  time- 
dependent  region,  we  must  solve  the  scalar  field  equation 


»  =  —  W=G  go3  *  )  -  0 

vCg  >B  *“ 


(3.3) 


G  =  det  g 


a3 


using  the  background  metric  (3.1)  for  that  region  [59].  From  the 
spherical  symmetry  of  $  ,  (3.3)  then  becomes  simply 


$ 

xy 


a 

2r 


($  +  $  )  =  0 
x  y 


(3.4) 


where  the  subscripts  here  denote  partial  differentiation. 

It  now  remains  to  find  the  initial  values  on  the  characteristics 
y  =  0  ,  x  =  xq  .  This  we  will  do  by  finding  exact  solutions  of  (3.3)  in 
both  the  static  interior  and  static  exterior  regions  (Region  I  and  Region 
II  of  Figure  2,  respectively),  and  using  the  fact  that  these  solutions  (but 
not  their  first  derivatives)  must  match  -  not  only  with  each  other  across 
their  common  boundary  (r  =  Rq)  ,  but  also  with  the  non-static  solution 
across  the  characteristics  which  bound  the  time-dependent  region. 


First  we  solve  (3.3)  in  the  static  interior  region.  In  this 
case  the  background  metric  is  flat  [60],  and  (3.3)  becomes,  in  the  case  of 
spherical  symmetry. 


r  $  +  2  $  =0  , 

rr  r 


which  has  the  general  solution 
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$  =  c1  4-  c2/r  .  (3.5) 

In  the  static  exterior  region,  we  have  the  Schwarzschild 
background  (1.1),  and  for  static  spherical  symmetry,  with  2m  =  1  , 

(3.3)  becomes 


(r-l)  $rr  +  2(1  -  i)  =  0  . 

The  solution  in  this  case  [61]  is 

$  =  c0  +  c,  tn  (— -) 

3  4  r 

The  arbitrary  constant  in  $  makes  no  difference  to  the  field,  so  we 
choose  c^  =  0  for  simplicity,  and  hence 

$  =  c^  In  (r-^-)  •  (3.6) 

At  this  point  it  is  convenient  to  assume  the  initial  radius 
of  the  shell  is  very  large  (R  >>  1)  .  Then  throughout  the  static 
exterior  region  (r  >_  Rq  >>  1)  we  have 

£n(— )  ~  .  (3.7) 

r  r 

In  particular,  in  order  that  (3.5)  and  (3.6)  match  up  across  r  =  Rq  >>  1  , 
it  follows  from  (3.7)  that 


c,  +  c0/R 
1  2  o 


/R 
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This  tells  us  that 
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(3.8) 


From  (3.8)  our  static  interior  solution  is  therefore 


$  =  -  , 
r 


(3.9) 


and  from  (3.7)  and  (3.8)  the  static  exterior  solution  is 


$  =  -  c  In  (^)  ~  - 
r  r 


(3.10) 


The  value  of  $  on  the  characteristic  y  =  0  (history  of  the  collapsing 
shell)  is  then  determined  by  the  solution  (3.9),  and  so  (choosing  c  =  1 
for  simplicity)  we  have 


4>(x,y=0) 


1 

r(x,y=0) 


(3.11) 


Similarly,  for  the  characteristic  x  =  x  =  2(R  -1)  +  2  £n(R  -1)  >>  1 

o  o  o 

(history  of  the  shock  front)  we  use  (3.10)  to  conclude  that 


$  (x=°°  ,y) 


1 

r  (x=°°,y) 


(3.12) 


If  we  define  \p  =  r  $  ,  the  characteristic  initial  value  problem 
given  by  (3.4),  (3.11)  and  (3.12)  takes  the  simple  form 


a 

4r 


ip(xt0) 


5 


(3.13) 


(3.14) 


' 
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§3.3  Magnetic  Dipole  Model 

We  now  turn  to  our  second  model  in  which  a  static  axi-symmet ric 

magnetic  dipole  of  moment  y  is  placed  at  the  centre  of  the  shell.  As 

2 

before,  it  is  assumed  that  y  <<  1  ,  i.e.  the  perturbing  field  is  weak 
enough  that  it  in  no  way  affects  the  geometry  of  the  system. 

Our  first  step  is  to  write  down  the  source-free  electromagnetic 
field  equations 


F^l  =  0  , 

v  * 


F  r  -I  =  0 

[yv,a] 


where  the  electromagnetic  field  tensor  is  given  by 


(3.15) 

(3.16) 


F  =  A  -  A 
yv  v,y  y,v 


(3.17) 


and  A^  =  (Aq ,A^ ,A^)  is  the  associated  potential  4-vector.  We  need 
concern  ourselves  only  with  (3.15),  since  (3.16)  is  satisfied  identically 
by  virtue  of  (3.17). 

As  in  §3.2,  we  want  to  solve  these  field  equations  in  the  time- 
dependent  region,  given  initial  values  on  the  characteristics  y  =  0  , 
x  =  xq  .  If  we  write  out  the  equations  (3.15)  on  the  Schwarzschild  back¬ 
ground  (3.1)  and  use  (3.17),  together  with  the  facts  that  the  field  is  axi- 
symmetric  and  gauge  invariant,  we  obtain  the  following  results  (see  Appendix 


B): 
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y  =  3  : 


a 


3  2 

( — ±&±\  _  o 

N  „  -!  „  Q  >  O  U  » 


sin  0  3,01  _  2  vsin  0'  ,2 

2r 


(3.18) 


M  -  2  :  A_  -  f(x,0)  +  g(y,0)  ,  f  and  g  arbitrary, 


(3.19) 


y  =  1  : 


y  =  0  : 


4(~  An  ,  )  n  +  -- j-1  -  -  (sin  0  A  )  =  1  (sin  0  A  )  , 

a  0,1  ,0  sm  0  0,2  ,2  sm  0  2,0  ,2 

(3.20) 

r2  1  1 

4(1T  Al,o\l  +  Sin~e  (sln  0  A1,2),2  =  iLTe  (sin  9  A2,l),2  • 

(3.21) 


As  before,  in  order  to  obtain  suitable  initial  conditions,  we 
must  find  corresponding  solutions  of  the  field  equations  for  the  static 
interior  and  static  exterior  regions.  Because  these  solutions  are  static, 
it  is  convenient  to  use  the  standard  Schwarzschild  coordinates  xV  =  (t,r,0,<j>) 
which  are  related  to  our  original  xV  =  (x,y,0,t{>)  by  (3.2).  The  corres¬ 
ponding  components  of  the  potential  4-vector  are  then  given  by 


so  that 


A  =  2 (A  -A  ) 
o  1  o 


-  (A  +A. ) 
a  o  1 


(3.22) 


Consider  first  the  static  interior  region  where  space-time  is 
flat  (a=l)  .  From  the  Maxwell  equations  with  magnetic  field  H  and 
vanishing  electric  field,  we  can  write 
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div  H  =  0  , 

f++j  7 


from  which  it  follows  that 


(3.23) 


H  =  curl  A  =  grad  ft  ,  (3.24) 

where  A  =  (A^A^A^  is  the  3-vector  potential.  By  virtue  of  (3.23) 
and  (3.24)  we  have 

V2ft  =  0  , 


which,  for  a  magnetic  dipole  of  moment  p  ,  yields  a  solution 


ft  =  - 


U  cos  0 


(3.25) 


Taking  contravariant  components  of  (3.24)  gives 


2  ipq  A  ij  n 

8  q,P  8  ,j  ’ 


(3.26) 


where  g  =  det  g  and  e1^5^  is  the  standard  permutation  symbol.  Using 

(3.25),  (3.26)  has,  as  one  of  its  solutions 


A^  =  A2  =  0  ,  A^  =  y  sin20/r  . 


(3.27) 


(In  addition,  A^  can  be  chosen  zero  since  there  is  no  electric  field). 
Returning  to  our  original  coordinates  xV  ,  a  solution  of  the  field 
equations  in  the  static  interior  region  is  now  found,  from  (3.22),  and 


(3.27),  to  be 
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y 


•  20 
sin  6 

r 


(3.28) 


One  can  easily  check  that  this  solution  satisfies  (3.18)  -  (3.21)  with 
a=l  ,  as  required. 


Turning  to  the  static  exterior  region,  we  must  find  the  corres¬ 
ponding  static  solution.  It  is  clear  that,  in  this  region  as  well,  we 
may  choose 


A 

o 


0 


(3.29) 


so  that  (3.19)  -  (3.21)  are  trivially  satisfied.  Let  us  define 


<Kx,y) 


•  2  Q 

sm  0 


(3.30) 


(Comparison  with  (3.28)  shows  that  ip  =  y/r  represents  the  appropriate 
static  interior  solution.)  From  (3.30),  we  can  reduce  (3.18)  to  get 


<1* 


a 


xy  2r2 


(3. 31) 


where  again  the  subscripts  denote  partial  differentiation. 


For  a  static  solution,  it  is  easily  shown  from  (3.2)  that 


,  a  .  ,  a  . 

ib  =  —7 -  ip  +  — 77  4> 

xy  4  rr  ,2  Tr 
J  4r 


(3.32) 


Combining  (3.31)  and  (3.32)  we  obtain 


x  2 
a  Ip  +  —77  Ip  =  —77  Ip  , 

rr  2  r  2 
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which  can  be  solved  to  give 


dr 

r3 (r-1) 


C  arbitrary  constant.  (3.33) 


If  we  assume  as  before  that  >>  1  ,  then  r  >>  1  throughout  the  static 

exterior  region,  and  (3.33)  becomes 


*  ~  ~  ,  (3.34) 

with  suitable  choice  for  C  .  This  result  is  therefore  identical  to  the 
static  interior  solution. 

It  remains  to  find  the  manner  in  which  these  two  solutions  join 

up  with  the  non-static  solution  across  the  characteristics  y  =  0  , 

x  =  xq  >>  1  .  As  it  turns  out,  the  jump  conditions  for  the  electromagnetic 

field  require  continuity  of  the  across  the  characteristic  surfaces 

(see  Appendix  C).  In  the  time-dependent  region,  then,  we  can  again  take 

A  =  A,  =  A_  =  0  ,  so  that  the  characteristic  initial  value  problem  is 
o  1  2 

given  simply  by 


<J>(x,0) 

«K°°,y) 


y 

r (x,0) 

y 

r(°°,y) 


(3,35) 


(3.36) 


Comparison  with  the  scalar  monopole  case,  (3.13)  -  (3.14),  shows 
that  these  two  collapse  models  are  strikingly  similar  [62]. 
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§3.4  Results  of  Numerical  Integration 

The  two  characteristic  initial  value  problems,  (3.13)  -  (3.14) 
and  (3.35)  -  (3.36),  can  now  be  integrated  numerically  over  the  time- 
dependent  region  -«>  <x<°°  ,  0  _<  y  <  °°  [63].  Some  results  of  these 

integrations  can  be  seen  in  Figure  3  (scalar  monopole  perturbations)  and 
Figure  4  (magnetic  dipole  perturbations). 

Several  general  remarks  apply  to  both  cases.  To  a  stationary 
external  observer  (r  =  const.)  the  field  appears  nearly  constant  for 
a  period  about  equal  to  the  Newtonian  free-fall  time  [64]  i.e.  down  to 
x  ~  0  .  The  epoch  x  ~  0  is  marked  by  the  fairly  sudden  onset  of  an 
oscillatory  decline  towards  zero.  On  the  horizon  r  =  1  itself,  the 
field  displays  a  similar  damped  oscillatory  behaviour  as  a  function  of  y  . 
A  free-falling  observer  close  to  the  collapsing  body  sees  no  decline  in 
the  field,  but  we  can  find  no  support  (at  least  in  these  idealized  models) 
for  a  suggestion  by  Ginzburg  [65]  based  on  a  quasi-static  analysis  for  a 
magnetic  dipole  perturbation,  that  the  field  becomes  infinitely  compressed 
against  the  body. 

With  the  knowledge  that,  in  both  models,  the  external  field 
decays  in  an  oscillatory  fashion,  it  is  natural  to  ask  what  becomes  of  the 
external  scalar  and  magnetic  field  energies.  We  can  answer  this  question 
by  making  use  of  a  certain  identity,  which  we  now  derive. 

From  (3.2)  it  is  immediate  that  both  the  scalar  and  magnetic 
equations,  (3.13)  and  (3.35),  are  of  the  same  form,  namely 
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^  as  a  function  of  retarded  time  (-x)  for  stationary  observers  with  radial 
coordinates  equal  to  1.57,  3,  and  5  Schwarzschild  radii.  The  scalar  potential  $  =  \\>/r  . 
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coordinates  equal  to  1.5,  2.5,  4,  and  6  Schwarzschild  radii.  The  vector  potential 
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ipxy  =  f  (x+y)  ip  .  (3.37) 

Using  (3.37),  then,  we  can  easily  obtain  the  identity 

4  %2  +  f*2)  E  ^  (4,x2  +  f,p2)  ■  (3-38) 

Physically,  (3.38)  is  just  the  explicit  form  of  the  law  of  conservation 
of  energy  [66]. 


Now  we  integrate  (3.38)  over  the  time-dependent  region 
-°°  <  x  <  °°,  0  y  <  °°  ,  and  obtain 


+  f*2) 


r°° 


dy  + 

x=°° 


(* 


X 


2 

fr) 


y=0 


dx 


(3.39) 
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f*2> 


dy  + 

x— — 00 


(t 


2 

fr) 


dx 


y=c 


The  left  side  of  (3.39)  can  be  integrated  for  both  the  scalar  and  magnetic 
cases  by  making  use  of  the  characteristic  initial  values  (3.14)  or  (3.36). 
Since  the  field  energy  is  conserved,  it  follows  that  the  two  terms  on  the 
right  side  of  (3.39)  represent  (up  to  a  constant  factor)  that  energy  which 
falls  in  through  r  =  1  ,  and  that  which  is  radiated  out  to  infinity, 
respectively. 


Consider  first  the  scalar  monopole  case,  for  which 

f  =  a/4r^  =  (r-l)/4r^  .  (3.40) 

If  we  integrate  the  left  side  of  (3.39),  using  (3.14),  (3.40)  and  the  fact 
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that  dx  -  dy  -  2  dr/a  from  (3.2),  we  obtain  the  value  1/4  .  In 
addition,  the  contribution  of  the  terms  fij/’  in  both  the  right  hand 
integrals  of  (3.39)  is  zero.  For  the  first  term,  this  is  true  because 
x  =  -oo  implies  r  =  1  ,  by  (3.2),  and  hence  f  =  0  ,  from  (3.40).  For 
the  second  term,  (3.40)  again  guarantees  that  f  vanishes,  since  r  =  00 
when  y  =  oo  [67],  Thus  (3.39)  becomes 


f°°  2 

dy  + 

2 

tr 

if 

Jo  y 

X=-oo 

a 

X 

— OO 

dx 


1 

4  * 


(3.41) 


The  two  remaining  integrals  in  (3.41)  can  be  evaluated  numerically,  with 
data  provided  by  the  computer  output  from  the  original  integrations.  The 
result  of  these  calculations  is 


OO 
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r  9 

if 

Jo  y 

dy  ~  .168  , 

X=-°o 

J 

i 

8 

-e- 

X 

t 

dx  ~  .080  . 

y=oo 

(3.42) 

Since  the  sum  of  these  integrals  is  very  close  to  ,  almost  all  of  the 
scalar  field  energy  (i.e.  99%)  is  accounted  for.  We  may  therefore  say  with 
confidence  that  about  two-thirds  (approximately  67%)  of  the  energy  falls  in 
through  the  horizon,  while  the  remainder  is  radiated  away  to  infinity. 

We  now  turn  to  the  magnetic  dipole  model.  This  time 

f  =  .  (3.43) 

2r  2r 

Using  (3.36)  and  (3.43),  the  same  approach  as  before  yields  the  value 

o  2 

3y  /8  for  the  left  side  of  (3.39).  Again  the  contribution  of  fip  to 

both  of  the  remaining  integrals  is  zero,  from  (3.43),  so  that  (3.39) 
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reduces  to 
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(3.44) 


Numerically  integrating  the  left  side  of  (3.44)  gives 


fCO 


'0 


X=-oo 


dy  ~  .360  y‘ 


00 


X 


dx  ~  .  015y‘ 

y=co 


(3.45) 


Here,  virtually  all  of  the  electromagnetic  field  energy  is  accounted  for, 
and  the  corresponding  proportions  which  pass  in  through  the  Schwarzschild 
radius  and  out  to  infinity  are  96%  and  4%,  respectively  [68]. 

§3.5  Discussion  and  Remarks 

In  the  previous  three  sections  we  have  considered  two  highly 
idealized  models  of  perturbed  collapse,  involving  weak  scalar  monopole 
and  magnetic  dipole  fields  anchored  to  thin  spherical  shells  collapsing 
with  the  speed  of  light.  We  have  seen  that  the  perturbing  field  decays  in 
an  oscillatory  manner,  with  most  of  the  field  energy  falling  into  the 
"black  hole"  formed  when  the  object  collapses  inside  r  =  1  .  A  summary 
of  these  results  appears  in  Table  1. 

We  conclude  this  section  with  several  remarks  worth  noting. 

(i)  In  the  electromagnetic  model,  only  4%  of  the  field  energy  escapes 
to  the  outside,  as  compared  with  33%  for  the  scalar  model.  The 
difference  here  is  presumably  due  to  the  fact  that  scalar  fields 
can  emit  monopole  radiation,  while  electromagnetic  fields  cannot. 
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TABLE  I:  Summary  of  Analysis  for  Two  Idealized  Models  of  Perturbed  Gravitational  Collapse 
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(ii)  Israel  [69]  has  described  the  oscillatory  nature  of  the  decay 
as  a  Ubathplug  effect”:  not  all  of  the  inward-falling  field 
energy  is  sucked  into  the  hole  at  once;  instead,  some  of  it  is 
turned  aside  by  field  pressure,  and  it  swirls  about  the  hole  on 
a  time-scale  several  times  larger  than  the  Schwarzschild  charac¬ 
teristic  time.  In  this  same  connection,  Price's  analysis  [70] 
establishes  the  existence  of  a  curvature-induced  potential 
barrier  at  r  ~  3m  which  prevents  the  escape  of  low-frequency 
waves  to  infinity.  Further,  the  potential  barrier  causes  the  low- 
frequency  waves  to  oscillate  -  once  for  monopole  fields  (t= 0) , 
twice  for  dipole  fields  (£=1),  etc.  [71]. 

(iii)  With  regard  to  the  second  (magnetic)  model ,  our  results  have  an 
important  bearing  on  the  question  of  observability  of  collapsed 
stars,  and  their  distinguishability  from  neutron  stars.  Unlike 
neutron  stars,  which  can  remain  intensely  magnetized  for  periods 
exceeding  10  years  [72],  a  collapsed  object  cannot  be  a  source 
of  synchrotron  or  other  radiation  requiring  the  agency  of  magnetic 
fields.  However,  accretion  of  interstellar  material  in  a  favour¬ 
able  environment,  e.g.  in  close  binary  systems,  could  produce  a 
strong  source  of  thermal  x-ray  bremsstrahlung. 

(iv)  Again,  concerning  the  magnetic  dipole  model,  a  glance  at  Figure  4 
reveals  a  curious  fact.  The  r  =  const.  curves  appear  to  have 
several  points  of  coincidence.  This  has  been  further  verified  by 
plotting  similar  curves  (not  shown)  for  several  additional  values 
of  r  ,  with  the  same  result.  We  have  found  no  explanation  for  this, 
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but  it  may  possibly  have  significance  regarding  a  resonant 
quality  of  space-time  [73]. 

(v)  A  similar  analysis  for  the  case  of  a  weak  gravitational  quadru- 
pole  perturbation  has  been  carried  out  by  de  la  Cruz  and  Israel 
[74].  The  results  lack  the  completeness  and  accuracy  which  our 
models  yield,  but  they  are  qualitatively  alike.  In  particular, 
only  a  small  amount  of  gravitational  energy  escapes  to  infinity. 

(vi)  Our  efforts  to  find  the  asymptotic  behaviour  of  the  external 
field  have  unfortunately  failed.  Price  [75],  however,  does 
succeed  in  showing  how  the  tails  of  the  radiated  waves  die  out 
for  large  t  .  He  considers  weak,  zero-mass  fields  of  arbitrary 
integer  spin  s  on  a  Schwarzschild  background,  and  shows  that,  as 
t  -*  °°  ,  a  multipole  moment  of  order  Z  decays  in  time  as 
Zn  t/t  (Z  >  s  >  0)  or  1/t  (Z  =  s  =  0)  [76], 
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CHAPTER  IV 

Event  Horizons  in  Static  Scalar-Vacuum 

Space-Times 


§4.1  Introduction 

Recent  interest  in  the  theory  of  gravitational  collapse  has 

raised  many  questions  regarding  the  existence  and  nature  of  event  horizons 

in  relativity.  Some  definite  results  are  known.  Israel  has  shown  [77] 

that  for  the  class  of  asymptotically  flat,  static  vacuum  fields,  only 

the  spherically  symmetric  Schwarzschild  solutions  with  m  >_  0  have  a 

regular  event  horizon  (r  =  2m),  and  [78]  that  for  the  corresponding 

1 

2 1  |  . 

electrovac  space-times , the  Reissner-Nordstrbm  solutions  with  m  ^  G  |e|/c 
are  the  only  ones  with  non-singular  horizons.  In  view  of  these  special 
cases,  it  is  therefore  natural  to  ask  whether,  for  arbitrary,  asymptotically 
flat  static  fields,  a  regular  event  horizon  is  destroyed  by  any  asymmetric 
perturbation  due  to  sources  within  the  surface  g  =  0. 

In  this  connection  there  has  been  some  recent  interest  in 
another  special  class  -  namely,  the  coupled  gravitational  and  massless 
scalar  fields  (where  by  "massless  scalar  field"  we  mean  a  scalar  field 
for  zero-mass  particles).  The  spherically  symmetric  solution  of  Janis, 
Newman  and  Winicour  (JNW)  [79]  has  the  interesting  property  that  the  event 
horizon  is  a  singular  point  in  the  space  no  matter  how  small  the  coupling 
constant  becomes.  Penney  [80]  has  suggested  that  this  surprising  result 
is  due  to  the  imposition  of  spherical  symmetry,  and  that,  by  considering 
asymmetric  solutions,  one  is  led  to  a  non— singular  horizon.  However,  his 
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example  in  support  of  this  contention  contains  an  error  [81],  and,  in 
fact,  his  solution  is  singular  at  the  horizon. 

In  this  paper  we  propose  to  clear  up  much  of  the  controversy 
about  event  horizons  associated  with  asymptotically  flat,  static,  massless 
scalar  fields  interacting  with  Einstein  fields.  Our  main  result  can  be 
stated  as  a  theorem:  every  zero-mass  scalar  field  which  is  gravitation¬ 
ally  coupled,  static  and  asymptotically  flat,  becomes  singular  at  a 
simply-connected  event  horizon.  This  theorem  immediately  obviates 
Penney’ s  search  for  a  non-singular  asymmetric  horizon. 

We  proceed  by  reformulating  the  given  conditions  in  terms  of 
the  geometry  of  the  surfaces  gQQ  =  const.  (§§4.2  -  4.3).  The  theorem 
is  stated  in  detail  (§4.4)  and  proved  (§4.5).  In  the  special  case  where 
we  neglect  the  gravitational  effect  of  the  scalar  energy  density,  the 
solutions  are  computed  explicitely  (§4.6).  In  §4.7,  we  discuss  some 
properties  of  the  singular  horizons. 

§4.2  Static  Fields 

This  section  deals  with  the  general  static  field.  We  want  to 
reformulate  the  Einstein  equations  as  conditions  on  the  geometry  of  the 
equipotential  surfaces  [82] . 

The  signature  of  the  metric  is  -4-H-.  Capitalized  Latin 
indices  run  from  0  to  3.  Three-dimensional  and  two-dimensional 
subtensors  are  distinguished  by  Greek  indices  (range  1-3)  and  by  lower¬ 
case  Latin  indices  (range  2-3).  Covariant  differentiation  with  respect 
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to  the  4-dimensional,  3-dimensional,  and  2-dimensional  metrics  is 
denoted  by  V,  a  stroke,  and  a  semicolon,  respectively. 


A  space-time  is  called  "static"  if  it  admits  a  regular, 

hypersurface-orthogonal.  Killing  vector  field  £,  which  is  time-like 
A 

(ZpK  <  0)  over  some  domain;  i.e.,  we  have 


V 

K 


A?B  +  V 


[AVC^B] 


0 


(Killing’s  equations), 
(hypersurface-orthogonality) . 


(4.1) 

(4.2) 


From  (4.1)  and  (4.2)  it  is  easy  to  show  that 

VV'2W  ■ 0  >  (4-3) 

where  we  define  V  by 

V  =  (-5acV'S  •  (4.4) 

Therefore,  throughout  a  simply-connected  domain  in  which  £  is  time- 
like,  (4.3)  enables  us  to  choose  t(x  )  such  that 

v'\  =  -3At  .  (4.5) 

This  allows  us  to  introduce  "static  coordinates"  x°  =  t,  xa,  where 
the  latter  are  any  three  independent  solutions  of 

£A8Axa  =  0  .  (4.6) 

A 

Using  (4.1),  (4.4),  (4.5)  and  (4.6),  it  is  straightforward  to  show  that, 
in  the  domain  where  <  0  ,  the  metric  can  be  put  in  the  form 
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a  1  (  1  2  3,,  a,  3  T72j  2 

ds  =  8ag  x  ,x  ,x  ^dx  dx  “  v  dt  » 


V  =  V(xb,x2,x2)  >  0 


(4.7) 


We  can  further  decompose  the  metric  form  (4.7)  if  we  define 


p  1  =  (v,  v|a)^2 

a 


(4.8) 


and  assume  p  d  vanishes  nowhere  in  the  domain  of  interest  [83] .  As 

intrinsic  coordinates  for  the  equipotential  2-spaces  V  =  const., 

1  2 

t  =  const.,  we  introduce  functions  0  ,  0  which  are  constant  along 

ot  8  ci 

the  orthogonal  trajectories  (g  3 ^0  3^V  =  0).  The  spatial  metric  then 
becomes 

ga^dxadx^  =  gab (V, 0)d0ad0b  +  [ p (V , 0 ) ] 2  dV2  .  (4.9) 


Suppose  n  is  the  unit  space-like  vector  normal  to  the 
equipotential  surfaces.  Then,  from  (4.8), 


n  =  p3  V  =  p  13xa(V,0)/3V  . 

a  a 


(4.10) 


In  addition,  we  let  e,  N  be  the  tangential  base  vectors  associated 

~(a) 

with  0  ,  so  that 


e“  v  =  3xa(V,0)/30a,  e^a)  =  gabe 
(a)  a 


(b)a 


=  30a/3xa 


(4.11) 


The  triad  (e,  .  ,n}  then  spans  the  3-space  at  each  point.  From  (4.9), 
~(a)  ~ 


(4.10)  and  (4.11),  or  by  simple  modification  of  (2.11),  it  is  immediate 


that 
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aft  ab  a  3  ,  a  3 

8  =  8  e(a)  e(b)  +  n  n 


(4.12) 


Using  (4.8),  it  is  easy  to  show 


2..  3  ~ 

p  Vi  Dn  =  -3  p 
a  3  a 


(4.13) 


The  extrinsic  curvature  of  the  2-space  V  =  const.,  considered  as  imbedded 
in  the  3-space  t  =  const.,  is  defined,  from  (2.4),  by 


*  a,  b  va  a 
6n  /S6  =  1^  e(a) 


(4.14) 


from  which 


T.  a  r  /JCQb  a  3  a  3  „ 

K  =  e. ,  N6n  /o0  =  e,  >.n  i0e,u>.  =  pe,  ^e,,  ^  Vi 

ab  (a)  a  (a)  a  3  (b)  (a)  (b)  a3 


(4.15) 


By  virtue  of  (4.12),  (4.13)  and  (4.15)  we  have 


\t  .(a)  (b) 

V  i  „  =  p  K  ,  e  e. 
a3  H  ab  a  3 


p~\p(eaC)n(3  +  egC>na) 


p  3(3p/3V)nang 


From  (4.10),  (4.11)  and  (4.14)  we  also  obtain 


(4.16) 


Kab  =  1  P_l38ab/3V  > 


(4.17) 


which  leads  to  the  related  formula 


sg’Vav  =  %2  pk  , 


(4.18) 


ab. 


where  g  is  the  2x2  determinant  of  g^b  and  K  -  g  Kab  is  twice 

From  (4.12)  and  (4.16)  it  follows  that 


the  mean  curvature. 


. 
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V'^  =  p  *K  -  p  3  3p/3V  . 


(4.19) 


The  Einstein  field  equations 


GAB  =  -8^Tab  > 


(4.20) 


(where  y  is  Newton’s  constant  of  gravitation  divided  by  c2) 
be  decomposed  [84],  yielding 


can  now 


J  =  8^T 


0  =  8ttyT 


ao 


-1 


G  .  =  -8ttyT  0  -  V  'L(V,  -  V'Y  g  )  . 

aB  a3  a3  y  sa3 


(4.21) 


The  relativistic  analogue  of  Poisson's  equation  is  then  given  by 


V1  V  =  47ryV(Ta  -  T°)  , 

y  a  o 


(4.22) 


which,  combined  with  (4.19),  gives 


p  28p/3V  =  K  -  4ttyVp  (Ta  -  T°) 

a  o 


(4.23) 


There  are  several  other  important  relations  which  we 
require  in  what  follows.  However,  their  derivations  are  somewhat 
tedious,  and  we  relegate  them  to  Appendix  D,  merely  listing  the  results 


here  for  later  use: 
1  1 


;b  1  „,b 


a  3  1  mA 


V_1g  23(g2VK^)/3V  =  -P  *“  -  -  8np(To6e"a)e"b)-  f  T^),  (4.24) 
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^(K  Kab-K2-R)  =  -8ttYT  Qnan6  +  p  1V  1K  , 
2  ab  aB 


(4.25) 


3  K  -  Kb  ,  =  -8ttYT  Qea,  +  p  2V  13  p  , 
a  a ;b  '  a3  (a)  K  a 


(4.26) 


Id  ABCD  ra3,  -2  -2  „ab  -4  -2  ;a  -6  -2.  ,,,,,2  n  . 

4^ABCDR  =  Ga3G  +p  V  KaKK  +2p  V  p.ap  +p  v  (9p/8v)  >  (4.27) 


ab 


3.b 

where  R  =  g  R  .  Equations  (4.17),  (4.23),  and  (4.24)  form  a  complete 
system  for  determining  the  evalution  of  g  ,  p ,  Kb  as  functions  of  V. 
Equations  (4.25)  and  (4.26)  are  involutive  constraints,  i.e.,  if  they 
are  satisfied  on  one  surface  V  =  const.,  they  must  be  satisfied 
identically.  The  last  equation  gives  the  (invariant)  square  of  the 
four-dimensional  Riemann  tensor. 


§4.3  Static  Massless  Scalar  Fields 


We  now  consider  a  static  scalar  field  $  for  zero-mass 
particles,  with  "scalar  gradient"  given  by 


a  a  =  -  (4tt)  3a$ 


(4.28) 


Since  $  is  static  we  have  a  =0.  The  scalar  equation  □$  =  0  can 

o 

be  written 


(_det  gAB}  2^a [  (-det  8AB)2§a3V]  =  °  * 


3 


(4.29) 


In  view  of  (4.7)  and  (4.9)  this  can  also  be  written  as 


V  13(Vg2i|/)/3V  =  -3a(pgVa)  , 


(4.30) 
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with  \p  defined  by 


a$/3v  =  pip  . 


(A. 31) 


For  a  massless  scalar  field  the  energy  tensor  is  given  by 


l  rn 

T  =  3  $9  $  -  —  e  2  9  $9  $ 

AB  A  B  2  8AB8  C  D  ’ 


(4.32) 


or  from  (4.28)  we  can  write 


4ttT°  =  -  \  a2  , 
o  2  * 


47TT°a  =  0  , 


,  rp0t3  a  8  1  a8  2 

4ttT  =  a  a  -  g  a  , 


(4.33) 


with 


a  =  -(4tt)^ (\pn  +e^a^9  $) 
a  ex  a  a 


2  _  a  .  2,  . ;aN 

a  -  a  a  =  4tt  (ib  +$  $  ) 

a  r  ;a 


(4.34) 


If  we  now  substitute  (4.33)  into  the  basic  equations  (4 . 23)- (4 . 26)  of  the 
previous  section,  we  are  led  to  the  following  complete  first-order 
system  for  determining  the  V-dependence  of  g  ,  $,  ip ,  p,  and  : 

ci  D  3. 


Geometrical  equation 


3gab/3V  =  2pKab 


(4.17) 


Static  Scalar  equations 


9$/9V  =  p\p  , 


(4.31) 
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1  1 


9  (Vg2^) /9V  =  -V3  (pgVa)  , 

cl 

(4.30) 

Gravitational  equations 


p  28p/3V  =  K  , 

(4.35) 

_  1  1 

V_1g  23(g2VKb)/3V  =  -p/b  -  ^  pR<$b  -  8ttyp$  .  <D;b  , 

cl  y  3.  A  d.  y  Si 

(4.36) 

Involutive  constraints 


^(KabKab-K2-R)  =  -47ry(^2-$>a$;a)  +  p_1V  ^  , 

(4.37) 

8  K  -  Kb  =  -8tty \b$  +  p  2V  bp 

a  a;b  ;a  ;a 

(4.38) 

The  following  result,  which  we  will  need  later,  is  obtained  by  contracting 
(4.36)  and  eliminating  R  by  means  of  (4.37): 


-1  .O  1  ?  0 

V8  (V  K)/3V  =  -p/  -  ^  pK  -pA  A  -  8t rypi/r  , 

y  Si  Z  SLu 

(4.39) 

where 


A  =  K  -  ^  g  ,  R 

ab  ab  2  ab 

(4.40) 

is  a  measure  of  deviation  from  spherical  symmetry. 


Finally,  we  combine  (4.27)  and  (4.35)  to  obtain 

IRABCDrABCD  '  GaBGa6+p"2v"2KabKab+2p"V2p;aP:a+p'V2K2  ' 

(4.41) 

ct  (3 

Although  an  exact  expression  for  the  quantity  G^G  could  be  evaluated 
from  (4.12),  (4.16),  (4.21),  (4.22),  (4.33),  and  (4.34),  it  is  sufficient 


•••• 
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to  note  here  that  all  terms  of  this  expression  will  be  positive,  with  one 
such  term  being 


9  9  9  .o  9 

48tt  y  (<T+$.  $’  ) 


(4.42) 


§4.4  Statement  of  Theorem 


In  a  static  space-time,  let  £  be  any  spatial  hypersurface 

A 

t  =  const.,  maximally  extended  consistent  with  £  £  <  0.  We  consider 

the  class  of  static  massless  scalar  fields  such  that  the  following 
conditions  are  satisfied  on  J  : 

(i)  £  is  a  "scalar-vacuum"  space  (i.e.,  free  of  matter  and 

sources  of  scalar  fields) . 

(ii)  J  is  regular,  non-compact  and  "asympotically  Euclidean". 

o i 

That  is,  there  exist  coordinates  x  in  terms  of  which  the 
metric  (4.7)  has  the  asymptotic  form 


8ag  =  6C6  +  0(r_1)’  Vas  =  °(r‘2)  ’ 

V  =  1  -  (m/r)  +  p  ,  m  =  const., 


n  =  0(r~2),  3  n  =  0(r  3),  3  3  n  =  0(r  4) 

0t  Oi  p 


>  (r  00 )  (4 .43) 


,  /x  a  3  >.2 

where  r  =  (ox  x  )  . 

a3 


(iii)  The  asymptotic  form  of  the  static  scalar  field  is 


$  =  (k/r)  +  C  >  k  =  const.. 


X,  =  0(r  2), 


3  c  =  0(r  3)  . 
a 


(r  -»■  00 )  (4.44) 
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(iv)  The  equipotential  surfaces  V  =  const.  >0,  t  =  const.,  are 
a  regular  family  of  simply-connected,  closed  2-spaces. 

(v)  If  the  greatest  lower  bound  of  V  on  £  is  zero,  then  the 
geometry  of  the  equipotential  surface  V  =  e  approaches  a 
limit  as  e  0  ,  corresponding  to  a  closed  regular  horizon 
of  non- infinite  area. 

(vi)  The  invariant  rAjBCDrABCD  is  bounded  on 


Theorem ;  There  is  no  non-trivial  static  space-time  which  satisfies 
conditions  (i)  -  (vi) . 


(Here  we  assume  y  >  0.  The  case  of  zero  coupling  is 
discussed  in  §4.6.) 


The  proof  of  the  theorem  is  presented  in  §4.5.  There  is  one 
trivial  case,  however,  which  can  be  quickly  disposed  of  here.  Suppose 
that  V  has  a  positive  lower  bound.  Then  the  maximally  extended  3-space 
£  is  complete.  Using  (4.29)  and  the  boundary  condition (4 . 44) ,  Green's 
theorem  implies  $  =  0.  Furthermore  (4.22)  reduces,  as  a  result  of 
comparing  (4.23)  and  (4.35),  to  Laplace's  equation  =  0  »  which, 

together  with  the  boundary  conditions  (4.43),  leads  to  V  5  1.  This 
means  space-time  is  flat,  and  the  theorem  is  established. 

We  may  assume  henceforth  that  V  comes  arbitrarily  close  to 
zero  on  £  .  The  equipotential  surface  V  =  0+  then  forms  an  inner 
boundary  of  £  .  Suppose  that  V  has  zero  gradient  at  some  interior 
point  P  of  J  .  Since  V  is  harmonic  (V^V  =0)  on  £  ,  P  would 
have  to  be  a  point  of  bifurcation  of  the  equipotential  surfaces  [85],  which 


. 
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contradicts  (iv).  Therefore,  by  (4.8),  p  ^  vanishes  nowhere  in  the 
domain  of  interest. 


We  conclude  this  section  by  recording  the  exterior  and  interior 
boundary  conditions  in  a  form  convenient  for  later  application.  For  the 
asymptotic  forms  (4.43)  and  (4.44),  we  find  from  (4.8),  (4.18),  and 
(4.31) 


,  2  -1 

r  p/r  m  ,  rk->2, 

2 

r$  k,  r  ip  -k,  as  V  1  . 


\ 


> 


(4.45) 


According  to  (vi)  and  (4.41),  the  regularity  of  the  manifold  at  the  inner 
boundary  V  =  0+  requires  that 

Kab  =  0(PV)>  P.a  =  0(p2V)  as  V  ->  0+  .  (4.46) 

It  follows  that  p  ^  is  constant  on  the  event  horizon: 


-1  12 

p  (0,6  ,0  )  =  1/pq  =  const. 


(4.47) 


ABCD 

In  addition,  since  the  curvature  scalar  R.^^R  is  bounded  every- 

ABCD  J 

where  on  £  by  (vi) ,  it  follows  that  the  expression  (4.42)  cannot 
become  infinite  anywhere  on  \  .  Hence  we  have  the  result  that  ip  and 
are  bounded  functions  throughout  J  . 


§4.5  Proof  of  Theorem 

In  this  section  we  employ  mainly  the  scalar  equations,  (4.30) 
and  (4.31), to  derive  integral  relations  which  enable  us  to  show  that 
the  trivial  solution  $  =  0  is  the  only  one  compatible  with  conditions 


(i)  -  (vi). 
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Let  F(V,$),  G(V,$)  be  (for  the  moment,  arbitrary)  differ¬ 
entiable  functions.  From  (4.30),  (4.31),  (4.35),  and  (4.18)  we  easily 

obtain  the  identity 
_1  1 

2  2—1  ? 

g  (9/9V)  [g  (VF(V,$)i|H-p  G(V,$)}]  =  A(V,$)p(f +$ .  $;a)  +  B(V,<W 

> a 

_  1  1  (4.48) 

+  p-19G/9V  -  Vg  2(Fpg2$;a)>  , 

> a 

where 


A  =  V9F/9$,  B  =  V9F/9V+9G/9$  .  (4.49) 

In  order  to  obtain  integral  conservation  laws  from  (4.48)  and  (4.49),  we 
require  that 


A  =  B  =  9G/9V  =  0  . 


(4.50) 


The  general  solution  of  this  linear  system  of  differential  equations  for 
F,  G  is  a  linear  combination  of  the  two  particular  solutions 


F  =  1,  G  =  1  , 

F  =  Jin  V,  G  =  -$ 


(4.51) 


Taking  these  values  in  turn  we  integrate  (4.48)  over  £  ,  i.e.,  we  form 

1 

2  i  2 

jj  j  (4.48)  g  dVd0  d6  ,  noting  that  the  integral  of  the  last  term, 
E 

being  a  2-divergence,  vanishes  when  taken  over  any  closed  2-space 
V  =  const.  The  results  express  the  equality  of  the  surface  integrals  of 
the  expression  in  square  brackets  above  over  any  two  equipotential 


surfaces  V  =  const.: 


. 
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/(ViJ>+p  1)ds  =  Cx 


(4.52) 


/[  (V  *n  V)i|i-p  HjdS  =  C2  , 


(4.53) 


2  12 

where  we  have  defined  the  element  of  area  by  dS  =  g  d0  d9 


As  an  immediate  consequence  of  (4.18)  and  (4.35)  we  have 


aCP  V)/9v  =  o  , 


(4.54) 


-1  A 

hence  if  we  form  Jjj  (4.54)  dVd6  <10  ,  we  obtain 


/  p  1dS  =  C,  . 


(4.55) 


Comparison  of  (4.52)  and  (4.55)  show  that 


/  Vt|>dS  -  C  -  C  =  C  . 


(4.56) 


We  can  now  evaluate  the  constants  C^,  and  by  integrating 

(4.53),  (4.55)  and  (4.56)  over  the  upper  boundary  V  =  1,  with  the 
help  of  the  boundary  conditions  (4.45).  We  thus  find,  as  integral 
conditions  on  the  lower  boundary  V  =  0+  : 


j  p  ^dS  =  47rm  , 


(4.57) 


V=0 


/  VipdS  =  -4irk  , 


(4.58) 


V=0 


+ 


/  [(V  Jin  V)ip-p  1$]dS  =  0 

v=o . 


(4.59) 


' 
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In  view  of  (4.47)  we  can  write  (4.57)  as 


S  /p  =  4Trm  , 
o  o 


(4.60) 


where  S  is  the  area  of  the  horizon  V  =  0. .  Since  S  is  non- 
o  +  o 

infinite  by  (v) ,  (4.60)  implies  that  is  also  non-infinite. 

We  now  consider  the  relation  (4.58).  Since  ip  and  the  surface 
area  are  bounded  on  the  horizon  V  =  0  ,  it  follows  that  k  =  0,  i.e.. 


/  VipdS  =  0  . 

v=o 


(4.61) 


+ 


In  addition,  the  boundedness  of  pip  =  3$/3V  and  $  ^  throughout  £ 
guarantees  that  $  itself  is  bounded  on  the  horizon.  This  fact, 
together  with  (4.61)  leads  to 


/  V$ipdS  =  0  . 

v=o . 


(4.62) 


We  now  return  to  the  identity  (4 . 48)- (4 . 49)  and  this  time  require 


A  =  V,  B  =  3G/3V  =  0  . 


(4.63) 


The  resulting  linear  differential  equations  have  the  particular  solution 


F  =  $  +  £n  V 


G  =  -$ 


(4.64) 


With  F,  G  given  by  (4.64),  we  thus  have  the  identity 

-II  .II 

g2 (3/3V)  [g2(VF^+p_1G)  ]  =  pV(i|;2+$>  $,a)  -  Vg2 (Fpg2$> ,a)  .  (4.65) 

)3.  5  d 


Integrating  over  \  (again  the  last  term  doesn't  contribute)  we  deduce 


' 
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the  inequality 


/  (VFiJH-p  1G)dS  >_  j  (VF^+p  1G)dS  .  (4.66) 

V=1  v=o+ 

From  (4.65)  it  is  clear  that  equality  in  (4.66)  holds  if  and  only  if 

►  (4.67) 

everywhere  on  £  .  Now  both  surface  integrals  in  (4.66)  can  actually  be 
evaluated  using  (4.64).  The  left  side  yields  the  value  zero  when 
integrated  with  the  help  of  (4.45),  as  does  the  right  side,  in  view  of 
(4.59)  and  (4.62).  We  conclude  that  (4.67)  must  hold,  hence  $  =  const, 
throughout  £  .  The  fact  that  $  vanishes  on  the  outer  boundary  V  =  1 
ensures  that  $  =  0  and  this  completes  the  proof.  It  is  worth  noting 
that  the  theorem  holds  even  for  a  regular  point  or  line  horizon  (S  =  0) . 

§4.6  Zero  Coupling 

We  consider  in  this  section  the  case  where  the  gravitational 
coupling  of  the  scalar  energy  density  is  neglected.  Our  problem  is  to 
obtain  solutions  of  the  vacuum  equation  G^^  =  ®  and  the  scalar  equation 
□$  =  0,  which  satisfy  conditions  (i)  -  (vi)  of  §4.4.  It  is  already 
known  [86]  that  the  only  vacuum  space-times  compatible  with  (i)  -  (vi) 
are  the  Schwarzschild  solutions  with  m  _>  0.  The  problem  thus  reduces 
to  finding  well-behaved  static  scalar  fields  defined  on  the  Schwarzschild 


il>  =  0 


$ 


a  =  0 


background : 
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8agdxadxg  =  (l-2m/r)  1dr2  +  r2(d02+sin2  ed*2) 


V  =  (l-2m/r)  . 


The  static  scalar  equation,  □$  =  0  ,  reduces  to 


Y2  3  ,v2  2  3<K  v2  3  /  .  n  3$x 

2  3r  (V  r  +  ~~2  ;  ”  39  (sin  6  30> 

r  r  sm  0 


2  2 

+  — J. _ 1_*  _  0 

2  2  2  u 

r  sin  0  3<f> 


.  (4.68) 

> 


(4.69) 


Separable  solutions  which  are  regular  on  the  axis  have  the  form 


$  =  R(r)P^(cos  0)e^^ 
n 


(4.70) 


where  R  satisfies 


2,  d2R 
(1-x  ) 

dx 


HR 

-  2x  —  +  n(n+l)R  =  0 


x  - - 1 

m 


(4.71) 

(4.72) 


For  general  n,  (4.71)  has  the  linearly  independent  solutions 

R=C.iP  (x),  R=C,  D  <x)  .  (4.73) 

in  z  n 

Of  these  solutions,  the  latter  are  unacceptable  for  all  n,  since  0  (x) 

is  singular  at  x  =  1  (i.e.,  at  the  event  horizon  r  =  2m).  On  the 

other  hand,  the  former  solutions  have  the  wrong  behaviour  at  infinity, 
except  when  n  =  0.  Hence  the  only  well-behaved  solution  satisfying 
(i)  -  (vi)  occurs  for  n  =  M  =  0  and  is  given  by 


' 
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*  =  cx  =  0  , 

from  the  boundary  condition  (4.44).  We  conclude  that  the  trivial  solution 
$  =  0  is  the  only  static  scalar  field  on  a  Schwarzschild  background  which 
is  well-vehaved  for  2m  <  r  <  00  . 

§4.7  Singular  Point  Horizons 

It  is  clear  from  the  previous  theorem  that  a  non-trivial 
static  scalar-vacuum  field  must  have  a  singular  event  horizon  (be  it  a 
point  or  otherwise),  and  thus  Penney’s  search  [87]  for  an  asymmetric 
solution  with  a  regular  event  horizon  is  unnecessary.  Our  theorem  also 
generalizes  the  work  of  JNW  [88]  who  show  that  all  spherically  symmetric 
solutions  always  have  singular  horizons. 

In  the  case  of  spherical  symmetry,  equations  (4.30)  and  (4.39) 
become  simply 

_1 

3(Vg2^)/8V  =  0  ,  (4.74) 

V8 (V_1K) / 9V  =  - -|  pK2  -  8uyp^2  .  (4.75) 

Now  (4.74),  (4.54)  and  (4.31)  can  be  solved  explicitly  with  the  help 
of  the  asymptotic  forms  (4.45)  to  give 


$  =  -km  £n  V  , 


(4.76) 


'<  J  I 


worfs  orfw  [88]  W  X  5o  'AiO'i  srfa  JisilaTaasg 


- 
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which  is  a  function  of  the  two  parameters  k  and  m.  Consideration 
of  (4.35),  (4.75)  and  (4.76)  enables  us  to  solve  for  p  and  K  as 
functions  of  V,  from  which  it  can  be  determined  that  p  -*  0  as  V  -*  0+. 
From  (4.57)  it  then  follows  that  the  area  of  the  inner  boundary  V  =  0+ 
is  zero,  i.e.,  the  horizon  is  a  point  (the  non-regularity  of  the  event 
horizon  means  that  (4.47)  is  no  longer  true,  but  the  integral  condition 
(4.57)  still  holds).  This  solution  is  precisely  the  JNW  solution  [89], 
which  has  a  singular  point  horizon  regardless  of  how  small  the  coupling 
constant  becomes. 

The  fact  that  the  horizon  is  always  a  point  in  the  case  of 
spherical  symmetry  leads  us  to  ask  whether  the  same  is  true  in  general. 
Penney ’s  example  [90],  although  not  regular  as  he  had  thought  [91], 
serves  to  show  that  there  are  asymmetric  solutions  whose  horizons, 
while  singular,  are  not  point-like.  His  axially  symmetric  solution 
has  the  line  element 


ds 


2 


e 


2v 


with 


v  =  —  2TTya 


(4.78) 


It  is  easy  to  check  that  the  horizon  (R  =  2m)  is  not  point-like.  This 
example,  along  with  (4.57)  suffices  to  show  that  the  function  p  need 
not  vanish  everywhere  on  the  horizon  V  =  0+.  We  can,  however,  show 


. 


* 
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for  a  general  scalar  solution  with  event  horizon  of  bounded  area,  that 
p  always  vanishes  at  least  locally  there.  From  (4.35)  and  (4.54)  it 
is  easily  shown  that 


/  P  2  dS  =  -/  p  1K  dS  . 


(4.79) 


From  (4.39)  and  (4.54)  we  have 


3 (p  1g2K)/3V 


-1  -1  2 
P  V  g  K 


1  1 

2  2  12  2 
8-ryg  iff  -  2  g  K 


-  g2I (£n  p)  ’a+P  2P 
» a 


(4.80) 


;  a  ,  .  . ab 

p  +A  ,  A  ]  . 
t  ab 


1  2 

If  we  now  form  jff  (4.80)  dVdQ  d9  and  use  the  Schwarz  inequality 
on  the  second  term  of  the  right  side. 


/  i|)2dS  i|l/  IMS]2  , 


it  is  straightforward  to  arrive  at  the  inequality 


^/Kp1dS<|/Kp1dS-^1[/  *dS]2  , 


(4.81) 


which,  in  view  of  (4.56)  can  be  written 


^  /  K  p"1  dS  <  |  /  K  p'1  dS  - 


8ttyCJ 

v2s 


(4.82) 


If  we  now  let 


. 

■ 
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X(V)  =  /  P-2  dS  , 


(4.83) 


then  (4.79)  and  (4.82)  combine  to  give 


d  r 1  dX,  C 

dV  dV^  —  8  *  C  -  const.  >  0 

VS 


(4.84) 


Since  S  S(V)  is  a  continuous  function,  infinite  only  on  the  outer 

boundary  V  =  1,  then  for  any  V1  <  1,  S(V)  must  attain  a  maximum  S.. 

i  M 

somewhere  in  the  interval  from  V  =  0+  to  V  =  In  this  interval 


we 


thus  have  1/S  >  1/S^,  hence 

—  M 


d_  ri  dX,  £]_ 

dV  dVJ  —  3  ’ 

V 


C'  =  C/Sr,  >  0 
M 


(4.85) 


Integrating  (4.85)  twice  we  arrive  at  the  inequality 


X(V)  >_  -C±  In  V  +  C2V  +  C3  , 


C1  =  C’/2  >  0,  C2  and  arbitrary  , 


(4.86) 


from  which 


X(0)  =  / 


-2 


dS  =  00 


V=0 


+ 


(4.87) 


The  integral  condition  (4.87)  along  with  (4.57)  shows  that  p  becomes 
zero  at  least  locally  on  the  inner  boundary  V  =  0+  (this  reaffirms 
that  for  spherical  symmetry  the  horizon  must  be  a  point).  It  is  also 
clear  from  (4.57)  that  the  region  of  the  horizon  on  which  p  vanishes 
must  be  of  zero  area.  In  other  words  the  gravitational  flux  p  =  | VV 


'i 


-  r  3,;  303 
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becomes  infinite  somewhere  on  the  horizon  because  some  finite-area  flux 
tube  shrinks  to  a  point  or  line  there. 


In  the  case  of  Penney' s  example,  comparison  of  (4.7)  and  (4.9) 
with  (4.77)  shows 


(4.88) 


2  2  2v  r  R  ^  ,2 

p  dV  ■  e  dR  • 


(4.89) 


From  (4.88)  and  (4.89)  we  immediately  deduce  that 


v  2  . 

p  =  e  R  /m 


(4.90) 


On  the  horizon  R  =  2m,  consideration  of  (4.78)  shows  that  p  vanishes 
in  the  equatorial  plane  0  =  tt/2. 


We  can  summarize  our  results  for  static  scalar  fields  as 

follows : 

(1)  Every  massless  scalar  field  which  is  gravitationally  coupled, 
static  and  asymptotically  flat,  becomes  singular  at  a 
simply- connected  event  horizon. 

(2)  In  the  case  of  spherical  symmetry,  the  singular  horizons 
are  points. 

(3)  If  we  allow  asymmetric  solutions,  it  is  possible  to  find 
examples  for  which  the  horizon  is  not  point-like. 

(4)  Assuming  the  horizon  has  bounded  surface  area,  the  gravitational 
flux  p  ^  always  becomes  singular  somewhere  on  the  horizon, 
due  to  the  shrinking  of  some  finite-area  flux  tube  to  a  point 


or  line. 


< 
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APPENDIX  A 

2 -Dimensional  Kinetic  Theory  of  a  Simple 

Relativistic  Gas 


Consider  a  2-dimensional  distribution  (for  example  a  spherical 
shell)  of  a  relativistic  ideal  gas.  Let  n(p)dp  be  the  number  of 
particles  with  momenta  (p,p+dp)  per  unit  area,  as  measured  in  the  rest 
frame  of  the  gas. 

The  baryon  density  of  the  gas  is  given  by 

00 

a  =  /  n(p)mAdp  ,  (A.l) 

o 

where  m  is  the  atomic  mass  of  a  constituent  particle.  The  total 

energy  of  an  atom  with  velocity  u  is 

_  1 

E  =  mA(l-u2)  2  ,  (A. 2) 

hence  the  total  surface  energy  density  is 

1 

00  00  —  — 

a  =  /  n(p)  E  dp  =  /  n(p)mA(l-u  )  dp  .  (A. 3) 

o  o 

The  pressure  exerted  by  the  gas  is  the  normal  component  of  the  mean 
rate  of  momentum  transfer  across  a  line  segment  of  unit  length,  and  is 

given  by 


P 


2  o 


n(p)pu  dp 


(A.  4) 


■ 


. 


80  - 


1  2 

where  the  factor  arises  from  averaging  cos  0  over  all  directions 
in  a  plane  and  the  momentum  p  of  an  atom  is  simply 

1_ 

p  =  m^u(l-u^)  ^  (A. 5) 

We  assume  that  pressure  is  a  monotone  increasing  function  of 

surface  energy  a,  and  discuss  the  special  limiting  cases  of  the 

equation  of  state,  namely  the  non-relativistic  limit  (P  <<  oq) ,  and 

the  ultra-relativistic  limit  (P  >>  a  ) . 

o 

In  the  non-relativistic  case  we  have  u  <<  1,  hence 
2  ~ 

E  ~  m  (1  +  u  /2)  and  p  ~  mu.  The  integrand  in  (A. 3)  then  vanishes 
except  for  very  small  p,  and  we  have 

00 

a  ~  a  +  yr  /  n(p)m  u^dp  ~  a  +  P  .  (A. 6) 

O  Z  A  O 

O 


Substituting  (A. 6)  into  the  isentropic  condition 

d(a/a  )  +  P  d(l/a  )  =  0  , 

o  o 


(A.  7) 


it  is  found  that 


P 


2 


a  a 

o 


(A. 8) 


The  dependence  of  P  on  a,  since  P  <<  cj^,  is  immediate  from 
(A. 6)  and  (A. 8),  hence 


.  - 


. 


. 


' 
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For  the  ultra-relativistic  case,  u  ->  1,  thus  E  ~  p,  from 
(A. 2)  and  (A. 5).  This  time  the  main  contribution  to  the  integrand  in 
(A. 3)  is  from  very  large  p.  It  is  then  immediate  from  (A. 3)  and  (A. 4) 
that 

oo 

a  Z  f  n(p)p  dp  ~  2P  ,  (A. 10) 

o 

i.e.,  P  a  a  ,  the  fully  relativistic  analogue  of  (A. 9).  Putting 
(A. 10)  into  (A. 7)  leads  to 


P  a  a  '  ,  (A. 11) 

o 

corresponding  to  (A. 8).  Equations  (A. 8)  and  (A. 11)  are  important, 
since  they  show  that,  for  an  ultra-relativistic  gas  exerting  pressure, 
the  adiabatic  exponent  y  =  3/2,  and  that  otherwise  3/2  <  y  <_  2. 

We  conclude  by  showing  that  a  -  2P  <  gq  for  a  simple 

2-dimensional  relativistic  gas,  interacting  only  by  collisions. 

1 

Since  p  =  mAu(l-u2)  2,  it  follows  from  (A.l)  ,  (A. 3)  and  (A. 4) 

that 

1 

oo  — 

2  2 

a  -  2P  =  /  n(p)mA(l-u  )  dp  <  oq 
o 


4 
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where 


With  x 


APPENDIX  B 

Electromagnetic  Field  Equations  on  a 

Schwarzschild  Background 


Consider  the  Schwarzschild  background  metric  in  the  form 


2  2  2 
ds  =  adxdy  +  r  , 


a  =  1 - 


(B.l) 


4  _ 


x  =  (r-1)  +  £n(r-l)  -  t  , 
y  =  (r-1)  +  £n(r-l)  +  t 


=  (x,y,e,tf>),  we  can  then  write 


0  a/2  0  0 


a/ 2  0 

0  0 

0  0 


0 


0 

0 


2  2 

0  r  sin  0 


,  [gyV] 


[Sv] 


-1 


(B .  2) 


(B .  3) 


We  are  concerned  with  the  electromagnetic  field  equations 


,yv 


v 


-i-  (/-G  FyV)  =  0 
/=G 


(B .  4) 


on  the  background  (B.l),  where 


■ 
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F  =  A  -  A 
yv  v,y  y,v 

Assuming  the  field  is  axi-symmetric  (independent  of  x~*  =  <j> 

2 

evaluating  G  =  det  g  =  -  \  r4sin2e  from  (B.3),  we  can 
(B.4)  as 


ar  sin0 


- -  (f  r2sin0  Fyj)  .  =  0 

Z  ,  J 


(j  =  0,1,2)  , 


or  equivalently. 


/0t  2  .  „  uo  icr 

(2  r  Slne  8  8  Fp0> , j  "  0 


Putting  y  =  3  in  (B. 6)  we  obtain 


/  a  01  .  ,  ,  a  W  v  .a  22 

2sin0  8  F3l\o  +  2sin0  8  F3o\l  +  ^2sin0  8  F32^,2  “  C 


which  from  (B.3)  and  (B.5)  simplifies  to 


2  .  ,  ot  >1  .  % 

A„  _  +  — —  ( — ; — —  A0  „)  _  =  0 


sin0  “3,01  '  2r 2  sin0  3,2  ,2 


In  a  similar  fashion,  (B.6)  with  y  =  2  becomes 


A  +  A 

1,02  0,12 


0  . 


Now  we  will  use  the  fact  that  the  field  is  gauge  invariant, 
unaltered  by  a  transformation  of  the  type 


=  A  +  $ 
y  ,y 


(B.5) 

and 

write 


(B.6) 


(B .  7) 

(B .  8) 

i.e. ,  it  is 


.  ' 
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Consequently 


Al,0  +  A0,l  Al,0  +  Ao,i  +  2$,oi 


Now  the  wave  equation  $  =  h(x,y,0)  can  always  be  solved;  in 

particular  we  can  choose  our  $  such  that 


$ 


,01 


2(A1,0  +  A0,P 


* 


so  that 


Ll,0  +  A0,l  ° 


(B .  9) 


can  always  be  made  to  hold.  In  view  of  (B.9),  then,  (B.8)  becomes 
simply 


A2,01  °  ’ 


from  which 


A2  =  f(x,9)  +  g(y,0),  f  and  g  arbitrary, 


(B.10) 


p ii^ia  1  ly  we  put  y  =  1  and  y  —  0  in  (B.6)  and  obtain  similar 


relations ,  namely 


(—  Fnl)  -  +  -At  (sine  F  )  -  0  , 

a  01  ,0  sinG  02  , 2 


(B. 11) 


1 

(—  ,  +  t  r  (sin0Flo)  =  0  . 

^  a  10  ,1  sm0  12  ,2 


(B . 12) 


vl  W..  v«  ohJalai 


- 
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It  is  easily  shown  from  (B.2)  that 


(—  )  =  (—  )  =  -2-r-~3 
a  ,0  {  a  \l  a 


(B . 13) 


and  hence  (B.ll)  and  (B.12)  become 


2r-3 


-  F„„  + 


2r' 


+  F, 


a  01  a  01,0  02,2 


+  cote  fQ2  =  o  , 


(B.14) 


a  F10  +  a  F10,l  +  F12,2  +  COt0  F12  ~  0 


(B . 15) 


Now  from  (B.5)  and  the  gauge  condition  (B.9)  we  have 


F01  2A0,1  »  F10  ~  2A1,0  * 


(B . 16) 


Finally,  using  (B.5),  (B.13)  and  (B.16),  we  can  put  (B.14)  and  (B.15) 
in  the  form 

r2  1  1 

4(a  A0,l\o  +  sine  (sin0  ^,2^,2  "  sine(sin0  A2,0),2*  (B’17) 


4(— —  A,  n)  +  — F— ■  (sin0  A..  9)  0  =  ~~  Q(sin6  A0  ..  )  .  (B.18) 

a  1,0  ,1  sine  1,2  ,2  sine  2,1  ,2 


The  four  key  equations  to  be  solved  for  A  are  therefore  (B.7),  (B.10), 


v 


(B . 17)  and  (B.18). 


It  is  worth  noting  that  the  left  sides  of  (B.17)  and  (B.18) 
involve  only  A^  and  A^,  respectively.  This  makes  it  particularly 
easy  to  see  that  A^  =  A^  =  =  0  trivially  satisfies  (B.10),  (B.17) 

and  (B.18). 


* 


.  ... 
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APPENDIX  C 

Jump  Conditions  for  the  Electromagnetic  Field 


We  are  concerned  here  with  finding  the  jump  conditions  for 
the  electromagnetic  field  across  an  arbitrary  3-surface  £  .  Suppose 
we  choose  a  3-dimensional  "cylinder"  R3  of  "length"  2e  and  interior 
R4  »  which  cuts  £  as  shown  in  Figure  5.  Applying  Stokes'  Theorem  to 
an  arbitrary  vector  quantity  Vy  ,  we  have 


(C.l) 


where  dxy  is  choosen  normal  to  the  3-surface  . 

We  now  let  Vy  =  FyV  ,  where  FyV  is  the  electromagnetic 
field  tensor  and  nv  is  arbitrary.  Then  (C.l)  becomes 


(C.2) 


The  first  term  on  the  right  side  of  (C.2)  vanishes. by  virtue  of  the 
electromagnetic  field  equations  (3.15).  In  the  limit  as  e  -*  0  ,  the 
volume  R^  becomes  sufficiently  small  that  the  right  side  disappears 
entirely,  while  only  the  ends  of  the  cylinder  R^  contribute  to  the  left 
side.  Therefore,  since  pV  is  arbitrary;  (C.2)  becomes 


F 

pv 


[F  dTP]  =  0  , 

pv  r 


90b1tu3-£  rs  HB  3201'JB  D.i:j£  ]f‘,!fO”l'3Di) I®  9ri3 


. 


ij'  ni.  o  I  ?.|  s  1  L  i  .  7  ,  :-t  ■  in‘> 


rr  aonis  ,97o5sM9ffT 


FIGURE  5 


The  3-surface  £  being  cut  by  a  3-cylinde 
One  dimension  is  suppressed  in  the  diagram. 


labn  I  o-C  b  v%d  .1  ud  sniad  7 
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where  dxy  is  normal  to  £  •  The  required  jump  conditions  are  thus 
given  by 


normal  to 


(C.3) 


As  a  particular  example,  we  consider  the  idealized  collapse 
model  of  §3.3,  in  which  the  two  characteristic  surfaces  of  interest 
are  given  by  y  =  0  and  x  =  x^  >>  1  .  Because  the  collapse  occurs 
with  the  speed  of  light,  these  surfaces  are  null,  which  means  that  a 
vector  normal  to  either  surface  is  also  tangent  to  that  surface. 


We  first  consider  the  surface  £  ,  given  by  y  =  0  .  For  the 

static  interior  field,  the  equation  of  the  collapsing  shell  is 


y  =  r+t  =  0  , 


from  which  a  normal  (hence  tangent)  vector  is 


Ny  -  ^  =  (-1,1, 0,0)  , 


(C.4) 


where  xy  =  (t,r,0,<|>)  .  Using  the  static  interior  solution 


—  —  —  „  —  u  sin  0 

A  =  A  —  A  —  0  ,  A  —  , 

o  1  2.  Jr 


found  in  §3.3,  together  with  (C.4)  and 


+  _  + 

F  I  =  A  -  A  | 
pv 1  y,v  v,y' 


(C.5) 


we  obtain 


zurll  si i.  enoiajtbijoo  q.  •  ;f,  boi rupai  sr 


0  3  1  BIOT  on  ii  Tb  97  9(1 W 


.  I  o3  J  smiorr 


■ 


:  b/iB  0  “  y  ■  navxg  sib 


fi'J  s  J  t  }.  sj!- 


io^dsv  ( Jrragn  3  son  »rl)  X*  b  /toirfw  moi5 


' 


n±E3do  9W 
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M=0, 1 ,2 :  F  N 


v 


y=3: 


yv 


MV 


=  o  , 


2 

M  sin  0 


nvi;  =  -  ^ 

I  zx  r 


(C.  6) 


Turning  to  the  exterior  field,  the  equation  of  the  shell  is 


y  =  (r-1)  +  £n(r-l)  +  t  =  0  , 


and  hence 


M 

dr 


=  (-  -,1,0,0) 


a 


a 


(r-1) 

r 


(C .  7) 


Using  (C . 3)  ,  (C. 5)  ,  (C.6)  and  (C.7)  we  easily  obtain  the  jump 
conditions  across  y  =  0  ,  namely: 


M=0,1: 

(A  , 
0,1 

“  A  )  |  y  - 

1^°  Zx 

0 

M=2 : 

(  Aft  2 
0,2 

b 

(A1,2  A2,l)^1  ° 

>  (C.8) 

M=3 : 

(A„  1 

-  A  /a) |r 

■  00 
\<  ' 

'V 

n 

+  -3 

M  sin  0 

3,1 

3,o  Zi 

dr 

t-3  • 

The  second  characteristic  surface  ^  is  given  by  x  =  >>  1  . 
If  we  now  carry  out  the  similar  analysis  for  this  case,  and  use  the  static 
exterior  solution 


M  sin  9 
r 


> 


found  in  §3.3,  the  corresponding  jump  conditions  turn  out  to  be 


?8  - 


' 


aiejdo 


90  - 


y=0 , 1 : 

(An  1 
0,1 

-Al,o>ll2=°  > 

\ 

y=2: 

0,2 

‘  S2,o)(^l£2  +  (X1,2  "  X2,l>ll2  ’  0 

)  (C.9) 

y=3: 

<V 

,  dA  +  .  2 

,7  /  \l+  3  y  sin  0 

+  Ao  /«)  v  ~  ,  v  =  o 

3’°  2-2  dr  h  r2 

• 

i 

The  jump  conditions  (C.8)  and  (C.9)  across  the  characteristic 
surfaces  y  =  0  ,  x  =  x^  >>  1  ,  guarantee  a  solution  in  the  time-depend¬ 
ent  exterior  region  of  the  form 

Ao  =  A1  =  A2  =  0  ,  a3  =  \p(x,y)  sin20  , 

where  ip  ~  y/r  on  both  characteristic  surfaces.  This  is  easily  seen 
physically,  because  the  characteristic  initial-value  problem  of  Chapter 
III  must  have  a  unique  solution  for  F  ,  and  since  the  assumption, 

Ao  =  =  a2  =  0  ,  made  in  Chapter  III,  satisfies  all  required  conditions 

(the  jump  conditions  (C.8),  (C.9)  as  well  as  the  field  equations)  it  must 


be  the  solution. 


7  •  +  A) 

o,t 

-I  ,i:3±ni  oiJeiia"  DSiBriD  sr.  1  <j;  moact  ^yilBDlsyriq 

■ 


■  nM  I  -'I 
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APPENDIX  D 
Imbedding  Relations 


We  first  derive  several  relations  for  imbedding  an  arbitrary 
hypersurface  £  in  an  (n+1)  -  dimensional  Riemannian  space.  Greek 
indices  run  from  0  to  n  ;  Latin  indices,  which  distinguish  quantities 
defined  on  the  hypersurface,  have  the  range  1  to  n  . 

Let  the  equations 

ot  ot  /  _1  . n  TTv  TT  /t\  i  \ 

x  =  x  (0  ,,,#,0  ;V)  ,  V  =  const.  (D.l) 


represent  an  orientable  hypersurface  £  with  tangent  base  vectors 


e^.^  =  3x  /30  , 


(D.2) 


and  unit  normal  n  : 


n*e,  . N  =  0  ,  n*n  =  e  (n)  =  { 


+1  n  spacelike, 


-1  n  timelike. 


(D.3) 


As  in  Chapter  II,  the  extrinsic  curvature  tensor  K  .  ,  and  the  intrinsic 


affine  connection  T  ^  are  defined  by 


e(a) 


e/  > ‘Se,.  ,/36  -  r  i  » 

~(a)  ~(b)  a,bc 


(D.4) 


'1 

,  >9qj  bni  niJBj  ;  <  03  0  wo  rt  «n  esoibnJt 
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from  which  we  obtain  the  Gauss— Weingarten  equations 


3e 


(a) 


-e(n)  K 


ab 


n 


+  r 


ab 


'(c) 


(D.5) 


(Note  that  in  Chapter  II,  £  was  a  timelike  hypersurface,  hence 

e(n)  =  +1)*  Then  (D.4)  and  (D.5)  ,  together  with  the  Ricci  commuta¬ 
tion  relations. 


—  — )e  )U 

39C  36b 


~  3  ~  Y 

^y  a  3x  3x 

o6y  6(a)  3  0b  3 9C 


lead  to  the  Gauss-Codazzi  equations 


(D.6) 


R  -  e .  “e,,  .‘e,  Ye,,*  =  R  ,  ,  +  e(n)(K,  K  ,  -  K  KJU)  ,  (D.  7) 

a3Y<S  (a)  (b)  (c)  (d)  abed  ~  dc  ad  ac  db 


a  3  Y  <5  _ 

a3Y<5  n  6(b)  6(c)  6(d)  *H>c;d  ^bd;c 


(D.8) 


Given  a  range  of  values  for  V  ,  (D.l)  then  gives  rise  to  a 

regular  family  of  hypersurfaces,  parametrized  so  that  01  are  constant 

ot  i  ot 

along  the  orthogonal  trajectories,  i.e.  n  30  /3x  =  0  .  It  then  follows 

that 

301/3x“  =  8lje(j)a  =  e(l)a  ,  (D.9) 

3xa(6,V)/3V  =  pna  ,  (D.10) 


where  p  is  defined  by 


' 


^cf  b  n  /■%  3l  q  9i9riw 
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n  =  e(n)p9  V  , 
a  ~  a 


or  equivalently 


P_1  =  (e(n)gaBOaV)06V)] 


-1 


where  p  vanishes  nowhere  in  the  region  of  interest  (see  §4.6 
follows  from  (D.12)  that 


2tt  3 

p  V.  n  =  -9  p 
a3  a 


Using  (D.4)  and  (D.10)  we  find 


3®(i)/3V  -  9(pn)/961*'  =  (3^p)n  +  p  e 


(j)  * 


from  which 


3gab/3V=  0/3V)(e(a).e(b))  =  2p  Kab  , 


9n/9V  =  -  e(n) (9  p)e 


(a) 


From  (D.4)  we  also  have 


tr  ~  /^Qb  a  3  a  3t, 

K  =  e,  .  •3n/90  =  e.  .  n  i.e,  ,  =  pe,  .  e.,.  V.  „  , 

ab  ~(a)  ~  (a)  a  3  (b)  (a)  (b)  a3 


which,  together  with  (D.13)  and  the  completeness  relation 


a3  ab  a  3  ,  /  \  a  3 

8  =  8  e(a)  e(b)  +  £(!l)n  n  • 


(D.  11) 

(D.12) 

It 

(D.13) 

(D. 14) 

(D.15) 

(D.16) 

(D.17) 

(D.18) 


£('  - 


b  r*  >w  (01. a)  bi  .  (£.a)  snleU 


3  r  ,  A  ,  8  jo  ds 

,(5>3  +  (d)B  (r,)  J  *  * 
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yields 


V,  =  p"1  K  ,e(a)  e(b) 


-2 


a3 


(a) 


(a)  a. 


abe  <xe  B  -  £<a>P  (3ap)(e  anB  +  6  gn) 


-  «(n)p  3(3p/3V)  n  n  . 
~  a  3 


(D.19) 


ab 

Then  for  the  mean  curvature  K  =  g  K  ,  we  have 

ab 


K  =  Pg0t3v|a3  +  e(~)p  2  9P(V»0)/9V  • 


(D. 20) 


If  we  now  form 


( (—  9  _  9  _9  \  \  a 

3V  361  '  391  3V  &  ’ 


and  use  (D.4)  ,  (D.5)  ,  (D.6)  ,  (D.14)  and  (D.16)  ,  we  obtain 


p  R  e .  “nynVe,  0  =  e(n)p  4-  g  9 K  P/9V  +  p  K  K,  p 
ayv3  (a)  (b)  ~  ;ab  ap  b  ap  b 


(D*  21) 


By  suitably  contracting  (D.7)  ,  (D.8)  and  (D.21)  ,  and  using 

(D.18)  ,  we  arrive  at  the  decomposition  of  the  Ricci  tensor,  R  0  =  Ry  „  , 

a3  aBy 

and  the  associated  Einstein  tensor,  G.=R0--^gnR,  with  respect 

*  a3  a3  2  toa3 

to  the  basis  {e/#N,n}  : 

~(i)  ~ 


ot  3  /  \  ab„  .  Tr  T,ab  T,2 

2G  .  n  n  =  -  e(n)g  R  u  +  K  UK  -  K  , 
a3  ~  b  ab  ab 


(D. 22) 


V^b/  "  3bK  - 


(D. 23) 


a_  3  _  Ti  _L  _j_  \  1/  V  _L  2V  P 


R  _e /  x  e/u  P  -  R  .  +  p  p  K  +  e(n)K  K  +  e(n)p  g  9K*79V  .  (D.24) 
a3  (a)  (b)  ab  ;ab  ~  ab  ~  ap  b 


(8.G)  <  0  .0)  gni^ofi  'Jaoo  Jct/a  vS 


n«  .  }  sl&d  9rf3  03 


- 


q*> 
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Let  us  now  consider  the  case  where  n  =  3  .  In  order  to  have 


Latin 


agreement  with  §4.2,  we  change  our  notation  so  that  capitalized 
indices  range  from  0  to  3  ,  while  Greek  indices  run  1  to  3  .  For 
the  spacelike  hypersurface  t  =  const.  (e(n)  =  -1)  it  turns  out  from 
(D.15)  that  the  extrinsic  curvature  y  V  ^3g  /8t  vanishes,  hence 

(D.22)  -  (D.24)  become  simply 


A  B  a3r, 
2GABn  n  =  g  Ra8  ’ 


p  A  B  n 

R,  n  e ,  .  =0  , 

AB  (3) 


RABe(a)  6 (3)  Ra3  +  V  V | a3 


> 


/ 


(D. 25) 


From  (D.25)  and  the  Einstein  field  equations 


gab  =  -  8l"<  tab  ’ 


we  readily  derive  the  following  decomposition: 


4  ga3  R  Q  =  8tty  T°  , 
2s  a3  1  o 


o  =  8  Try  T  , 

ao 


\ 


)  (D. 26) 


Ga$  =  ‘  8,nr  Tag  '  V_1(V|aB  "  V  'Vc.fS5  ’ 


We  are  also  interested  in  the  2-space  V  -  const. ,  considered 
as  imbedded  in  the  3-space  t  =  const.  This  time  we  use  the  imbedding 
relations  (D.22)  -  (D.24)  ,  where  n  =  2  ,  e(n)  =  +1  ,  and  lower  case 


. 
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Latin  indices  have  the  range  2,3  . 


Let  X  ,  Ya  ,  and  Z  ^  be  the  right-hand  sides  of  (D.22)  , 
(D.23)  ,  and  (D.24)  ,  respectively.  Then,  by  virtue  of  (A. 17)  (i.e. 

from  (D. 4)  ,  (D.9)  ,  and  (D.10))  we  have  the  following  identities: 


g  13(gX)/3V  =  PZab(Kab  -  Kgab)  -  p  1(p2Ya).a 
1  1 

g"  2  3(g2Y  )/3V  =  [p(6b  Za  -  Zb )  ]  +  p_1(p2X) 

a  a  c  a  ;  b  ;  a 


>  (D. 27) 


where  g  is  the  determinant  of  the  2-space  metric  g  ^  .  Thus,  if  (4.17) 
and  (D.24)  are  considered  as  a  system  of  first  order  equations  for  find¬ 
ing  g^  and  as  functions  of  V  ,  then  (D.22)  and  (D.23)  are 

"involutive  constraints",  i.e.  if  they  are  satisfied  on  one  surface  V  = 
const.,  then  from  (D.27)  they  must  hold  identically. 


From  (D.19)  ,  (D.24)  and  (D.26)  ,  we  obtain 


V  1  g  2  3 (g2  V  K  b)/3V  =  -o  ;b  -  -  P  R  6  b 


(D. 28) 


-  8Tryp(Taj3e(a)  e(b)6  2  g^)  , 


ah 

where  R  =  g  R  ,  ,  while  (D.22)  and  (D.23)  ,  respectively,  combine 

ab 

with  (D.19)  and  (D.26)  to  give 


i(K  Kab  -  K2  -  R)  =  -8wy  T  n"nS  +  p  *V  h  , 
2  ab  ap 


(D. 29) 


3  K  -  R  b  ,  =  -8tty  T  e.  “n3  +  P  “V  13  p  . 
a  a  ;b  a6  (a)  a 


(D. 30) 


■ 

a  ’  ^io  ms 'ic.  ;a  £  tei.  .fci  isb  no:>  k.C 

3o  eno  rwi  pS  line  §n/ 

(Ti.fl)  .oil  nsril  «.:)Bno3 


. 
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Finally,  from  (D.7)  ,  (D.8)  and  (D.21)  ,  with  suitable 


changes  in  notation,  we  find 


\  R  RABCD  =  G  GaB  +  V  2V,  V'a3  , 
4  ABCD  a3  a3 


which  together  with  (D.19)  ,  yields 


1  ABCD  a3  -2-2  ab  0  -4-2  ;a 

T  RADn>R  =  G  G  +pV  KUK  +  2p  V  pp 
4  ABCD  a3  ab  ;a 


+  p  6V  2  (8p/3V)2  . 


(D 


.31) 
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denotes  covariant  differentiation  in  the  timelike  hypersurface, 
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9u  9 v  9v  9u  ~  0Ay  9u  9v 
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2  2 
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of  convenience.  The  reasons  for  this  will  be  made  clear  in 
§3.2. 
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2  _  2  ,  ,  2  ,  2^2 
ds  =  -dt  +  dr  +  r  dfi 

[61]  The  general  static  solution  of  □$  =  0  on  a  Schwarzschild  background 

is  found  in  §4.6. 


[62]  That  the  similarity  of  (3.13)  and  (3.35)  is  not  a  coincidence  is 

demonstrated  by  Price  (ref.  18).  He  shows  that,  in  fact,  similar 
forms  occur  for  weak,  zero-mass  fields  of  arbitrary  integer  spin 
on  a  Schwarzschild  background. 


[63]  These  numerical  integrations  were  carried  out  on  the  University  of 

Alberta  IBM  360/67  system,  by  Mr.  R.  Teshima.  In  order  to  obtain 
equally-spaced  data,  we  employed  the  transformation  x  =  tan  X  , 
y  =  tan  Y  ,  so  that  the  time-dependent  region  was  given  by 


7T  IT 

' 2 -X-2 


0  <  Y  <  f 


Over  that  region,  the  computer  produced  a  1000*1000  grid  of  values 
for  ip  and  r  . 
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o 
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[66]  In  other  words,  (3.38)  can  be  obtained  from  the  conservation  law 

1/2 

((-g)  T  £^)  ^  =  0  ,  where  is  the  time-like  Killing 

vector  (see  4.2)  and  T  is  the  scalar  or  electromagnetic 
energy  tensor. 
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wise  we  would  have  true  singularities  on  the  horizon  (contradicting 
that  the  singularity  there  is  merely  an  apparent  one)  or  at 
infinity  (in  which  case  the  field  would  not  vanish  asymptotically). 
See  also  Novikov  (ref.  16). 


2  2 

[68]  The  actual  total  radiated  energy  in  this  case  is  not  3y  /8  but  y  /4 


(i.e.  the  constant  factor  here  is  — ) .  We  thus  have 


9  2  r°°  [ 

I,  +  I0  =  y  /4  where  I -=-=-/  ip 

1  2  1  3  rx 


y=c 


dx  and 


t  -  2  r  2 
Z2  "  3  4)  \ 

I1  =  .010  y2 


X=-oo 


dy  .  The  numerical  results  are  then 


I  =  .240  y  .  The  difference  between  the  magneto- 


1  9 

static  energy  -i-  y  initially  present  outside  r  =  1  ,  and  the  total 

2 

radiated  energy  y  /4  ,  is  exactly  accounted  for  as  the  work  done 
by  the  difference  of  the  electromagnetic  stresses  on  the  two  sides 
of  the  shell,  and  goes  into  increasing  the  kinetic  energy  of  the 
collapsing  body. 
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